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A Variational Principle of Hydromechanies 


S. DRoBot & A. RYBARSKI 


Communicated by C. TRUESDELL 


I. Introduction 


1. The purpose of this paper is to reveal the part played in variational princi- 
ples of hydromechanics by a certain group of infinitesimal transformations of 
the fields of density and velocity. This group, called here the hydromechanical 
variation and originating from some elementary requirement concerning variation 
of matter, seems to be essential in variational principles of hydromechanics, not 
only in deducing from them the equations of motion, jump conditions included, 
but also in studying the background of the conservation laws, the vorticity theo- 
rems included. 

So as to obtain such conservation laws, here we shall deal only with extremal 
variational principles, which assert that the equations of motion render the action 
an extremum. (They do not include flows in which dissipative forces act, e.g., 
flows of a viscous fluid.) The action will be taken as an integral, over a given 
time-space region, of any given Lagrangean depending on the fields of density 
and velocity only, these being considered as functions of the time-space point. 
The Lagrangean will not be a priori restricted to the usual, form, vzz., the dif- 
ference between the kinetic and the potential energy. However, it is not our 
aim to find more hydromechanical cases which can be included by these princi- 
ples or to study their range of relevance. We desire rather to ascertain which 
mathematical facts cause the equations of motion to follow from the variational 
principle, independently of the particular form of Lagrangean assumed; which 
facts lead to the conservation laws; and — as will be shown to be connected — 
which ones particularize the form of the Lagrangean. 

Namely, by suitably embedding the extremal principle of hydromechanics 
into the general framework of the calculus of variations, it will be proved that 
the invariance of the action with respect to a Galilean subgroup of the hydro- 
mechanical variation leads to the conservation laws of matter, of impulse and 
of angular momentum and restricts the form of the Lagrangean; however, the 
energy equation and — as seems to be more interesting and new — the vorticity 
conservation laws depend neither on this Galilean subgroup nor on the particular 
form of the Lagrangean but are subject to another intrinsic subgroup of the 
hydromechanical variation. 

Of course, not all results of this paper are claimed to be new. In appropriate 
places, some of the numerous existing related works will be quoted, but the list 
of them is by no means complete. We are greatly indebted to Professor C. TRUES- 
DELL for valuable discussion of the first draft of this paper. 
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2. In order to introduce the basic idea of the hydromechanical variation, 
we start by recalling that when in variational principles of hydromechanics we 
use Eulerian or spatial variables, some peculiar features arise which are not so 
evident when using Lagrangean or material variables. As in many earlier works, 
we find it instructive to compare the variational principles for continuous media 
with the principles valid in the mechanics of discrete material points. 

As typical of extremal principles for a system of N discrete material points 
with masses m,, @=1, 2,..., N, coordinates x,, 7=1, 2, 3, and potential energy 
U, let us consider HAMILTON’S principle: 


ty 
as 6farf Sm (GAs (GBs (GAY) — ul} =0 
to ves 


for all variations 6x‘(¢) vanishing at the ends of the time interval 4,</<4. 
We could, indeed, consider a more general integrand depending arbitrarily on 
all the variables occurring, but this would not illumine the questions in which 
we are interested. 


In hydromechanics formulated exclusively in terms of the Lagrangean or 
material variables a, b,c, the principle is merely an evident analogy to (1.1)?. 
It suffices to take instead of the finite sum over the indices a the integral over 
a given region V, of the triplets a, b,c; instead of the masses m, the element 
0 (a, b,c) dV), in which @, denotes the density of the fluid at the time /); instead 
of the coordinates x! the functions x’ (é, a, 0, c)yt=4, 2, 3; and to require, with 
a suitably chosen potential energy U, that 


(1.2) ofa {ff ants ele | (ary eo ul = 


for all variations Ox'(t, a, b,c) vanishing at the ends fy, 4, of the time interval. 
Guided in the same manner by the more general form of (1.1) with any given 
Lagrangean, we could obtain a corresponding generalization of (4.2)?. 


It is a different matter when we use exclusively the Eulerian or spatial 
variables and regard the velocity v(¢, x1, x2, x8) and the density 9 (t, x1, x?, x3) 
as fields defined at the point (x1, x?, x8) in a given spatial region V; and at the 
time ¢ in a given interval fj<¢<¢,. Should we try, purely formally, to write 
down the continuous analogue to (1.1) by taking, e.g., the potential energy «¢ as 
a given function, and by requiring 


t : 
(1.3) 6 J dt fff dV.eov%— e) =0 

fy Vi 
for all variations dg and 6v vanishing for f), 4, and on the boundary S, of the 
region V,, we should get from such a “principle” v= 0 as equations of “motion”. 
This cannot be avoided by taking a more general form of the Lagrangean, provided 
that it depends on ¢, x1, x?, x3 and on the fields 0,v only. The reason, of course, 


is that not all fields 9,v should be allowed in competition. It would seem to be 


1 Cf. e.g. HELLINGER [1914], p. 656, formula (25). 
* Cf. HELLINGER [1914], p. 656, formula (26) 
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enough to take the continuity equation 
00 3 a 
(1.4) =; +div(ov) =0 


expressing the conservation of matter, as a side condition to (1.3), or, what is 
mathematically the same, by introducing some Lagrangean multiplier g (é, x1, x2, x3) 
as an additional function to be determined from the principle, to require that 


as) afar fff ar{(Ler—) +9 (28 + aivies))} =o 


for all variations de, dv and 6@ vanishing on S, and for fj, t;. However, the 
equations of motion which follow from such a principle, viz. 

ao a ec a o = 28 1 div(od) = 

7? ihe Ope (Vera F + div(ev) =0 
pertain only to irrotational flow, gm being the potential of the velocity?. This 
result would not be improved by taking instead of the primary term $ov?—<« 
any other function depending on ¢, x}, x?, x8 and 9, d only. 

To obtain more general equations of motion one introduces‘ three additional 

MIMCUONS. OU 0", 24) AL, We, Xie) Le, Meee) and requires that 


(1.6) ) 9 fae avi{(> ov —e) + (52 + div(o#)) +4(36 +i grad yi)! =o 


for all variations do, 6v and dq, 6A, du vanishing on S, and for fp, ¢,.. The six 
equations which follow from this principle are 


1 2 | Ge Oy OM =_— — a 
Boge et set Midas Sava 
an coins BA (oma) . oe 
a, + ¥egradu =0, ay + div(ov) =0. 


The first of them is BERNOULLI’S equation; the next two are CLEBSCH’S equa- 
tions, and A =const, 4 =const are equations of the vortex lines®. By eliminating 
the three functions 9, A, “, we get the usual equations of continuity and momen- 
tum for the primary fields 9,v. Plainly the principle (1.6) can be generalized 
by taking instead of the term $ov?—e any other function depending on #, x1, x?, x8 
and 0,0 ; however, the additional terms will remain the same. 

It seems worthwhile to discover the part played by these additional terms, 
the introduction of which has been justified heretofore rather by the final success 
in getting formally the correct equations of motion than by any interpretation. 


3. To this end, let us analyse another approach, in which, however, both 
material and spatial variables are used. Namely, the continuity equation is used 


3 Cf. BATEMAN [1930], pp. 816—825. 
4 CLEBSCH [1859], pp. 1—10. Cf. also BATEMAN [1944], pp. 164— 166. 
5 Cf., e.g., LAMB [1906], pp. 233—234. 
Oe 
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in a material form®, : het 
2 é t, a,0,C 
(1.7) Q(t, #, %, #) = Oo 5G at a, xi)’ 


so that 


[ff odmdxdx=fff odadbdc 
Yi YN 


where V, is the region of places (x1, x?, x*) occupied at time ¢ by the particles 
a, b, c constituting the region V; taking a suitable expression I for the potential 
energy, one states, e.g., HAMILTON’s principle by requiring that 

ml =o 


by 
dx1\2 dx? \2 dx? \2 
a 0 far ff far yel| a eae { 
to Y, 
for all variations 6x’, i=1, 2,3 which vanish at the ends of the time interval, 
provided, however, that the total variation of the mass satisfies the condition 


(1.9) Aff eaVi=o0, 

which leads? to the following condition for variation of the density: 
a OO 

(1.10) Ae= Crags ) 


where the summation convention is used. While spatial variables are used in 
ax 

ay 

, which implies use of material variables. The principle can be generalized 


(1.10), in evaluating the variations of the velocities in (1.8) one takes a 
d 6x 
dt 
by taking as integrand an arbitrary function of the same variables*, the condi- 

tions (1.9) and (1.10) being retained. 


Just these conditions (1.9) or (1.10) will be proved to be essential in deducing 
the equations of motion from the variational principle. In fact, a special use 
of them is implicitly included in the introduction of the functions A, m, into 
the principle (1.6). 


In this paper we shall deal with the fields 0,v exclusively, using only field 
variables ¢, x1, x, x3 and not introducing additional functions, and we shall 
impose some conditions on the variations of the fields which constitute a gener- 
alization of (1.10). The physical background of these conditions can be found 
out by analysing more carefully some assumptions adopted, usually tacitly, in 


principles appropriate to discrete material points. Namely, one always assumes 
that 


(44) 0b = 0. andyldm 0), 


° Cf. HELLINGER [1914], p. 609, formula (7); also LICHTENSTEIN [1929], pp. 361 to 
365, particularly formulae (7), (8), (10) and (15). 

’ Cf. LICHTENSTEIN [1929], p. 343, formula (78); also HELLINGER [1914], p. 609, 
formula (8’), where a different notation is used. 

* In gas dynamics the internal energy is assumed to depend not only on the 
density but also on a third variable, such as the temperature. The problem of formulat- 
ing a purely spatial variational principle for such gas flows is one of classical difficulty. 
C. TRUESDELL has informed us that this problem has been solved recently by C.C. Lin, 
correcting insufficient earlier formulations by Eckart and HERIvEL. 
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and these are counterparts of (1.10). Further, one assumes that 


dm 
AeA sare 
(1.12) Ao he 
to which the continuity equation (1.4) or (1.7) corresponds. These conditions 
are logically independent and, although elementary, by no means trivial. We 
shall employ them in generalized form. Instead of the pair of assumptions (1.11) 
we shall adopt the continuous counterpart to the following: 


(1.13) 6t=0O implies dm=0O, 


and instead of (1.12), the continuous counterpart to the condition 
(1.14) (SF) Ot 


which makes possible inclusion of variational principles for points with variable 
masses and, in hydromechanics, flows with sources. Formulation of counterparts 
to the conditions (1.13) and (1.14) in terms of field variables leads to the definition 
of the hydromechanical variation in § II,1 below. 


4. We now explain the fields we shall use. In the Euclidean-Galilean four- 
dimensional space X, of the points x with the coordinates x*, « =O, 1, 2,3, by 
taking x® as the time ¢ and x’, 1=1,2, 3, as the space-like coordinates, we 
consider the four-dimensional vector fields #(x) with the components #*(x), in 
which ° is the density o, and #’, 1 =1, 2, 3, are the impulses 9 v', where v' denotes 
the component of the velocity v. 

The notion of the mass contained at the time ¢ in the volume JV, will be 
extended as follows. Let o* be any 3-dimensional hypersurface contained in X,, 
and let dt, denote the oriented element on o”’, 7.¢., 


(4:45) dq = Capys de dx) dx® 


where ¢,8,5= *878 is Ricci’s symbol, having the value 0 if two of the indices 

are equal, +1 if the permutation of the indices is even, and — 1 if the permuta- 

tion is odd, and where dx?, dx", dx°, are three linearly independent vectors 
1 


lying on o*; thus the vector dt, is normal to the hypersurface o*. The mass con- 
tained on o* will be represented by the integral fdr, p%, called also the complete 


matter-flow. In the particular case when the hypersurface o* is the space-like 
three-dimensional volume V, we have dt,»=dV, dt;=0, 1=1, 2, 3, and 


(1.16 farp*= [ff eaV 


reduces to the usual mass. More generally, when o* is closed and consists of 
V,,, V, and of the moving two-dimensional boundary S, of V, for fy<t<4, then 


ty = 

(1.47) $drap*= fff oaVe— IM eaV+ fat ff of a5, 
co Viz to to Sy 

this identity being known® as the matter-balance for the moving region. In the 

general case, when the hypersurface o* may be open, the complete matter-flow 


® Cf. CARTAN* [1922], pp. 36—37/. 
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represents a generalization of the notion of the mass contained in o*. By using 
Gauss’ formula 


(1.18) 6dt,p° =J adr 0,P*, 
in which t denotes any four-dimensional region contained in X,, and 7° its 
boundary, 0,= 6/0 x%, we see that 
2 s 00 . > 
(1.19) 6, p= 3; + div le?) 
represents the density of the sources of matter existing in T. 


The action W will be assumed to have the form 
(1.20) W=fdc Llp), 


in which the Lagrangean L is any given function depending on x and #(x) only. 
This form includes, e.g., the usual one for barotropic flow of an inviscid fluid, 
viz., as the special case when 

(1.24) L = sap {(e")? + (67)? + (69)9} — 26) — PU(A), 

where ¢ is the potential energy, which depends only on #°=g, and U is the poten- 
tial of the external forces. 


II. The variational principle 


1. In this section we formulate a variational principle of hydromechanics 
and deduce from it the equations of motion, whatever the assumed form of the 
Lagrangean may be. In this way we hope to show which facts are responsible 
for the derivation of these equations only. 

First we prepare some notions and formulas concerning variations. Let ‘8 
be a functional space of vector-valued functions p(x) with the components *(x) 
supposed to be sufficiently regular in X,. Consider the following infinitesimal 
transformations of the spaces X, and ‘8 into themselves: 


(2A) Gime t= OW Ae 

and 

(2.2) p=p+0(x). 

More exactly, 

(2.3) Ox(x) =e&(x)+0(e), dp(x) =ex(x) + 0(e), 


where € and z are arbitrary functions belonging to the space $8, e is a scalar 
parameter, and o(e)/e+0 as e->0. The functions d(x) are called the local 
variations of the field ~, and the functions Af given by 


(2.4) Ap = bf + agp - dx8, 


where @,=0/0x*, are called the total variations of bp. Yo within infinitesimals 
of first order we have 


(2.5) Ap = B(x) — p(x). 
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The operation 


(2.6) A fenpaelt. pone at ‘ 


where o* denotes the hypersurface obtained from o* by the transformation (2.1), 
will be called the total variation of the complete matter-flow. We have the 
identity 


(27) A Prac p* Tie {5p* — a, (p" bx* — p* bx?) — ,p°- dx", 
which can be proved by calculating the right-hand side of (2.6): 


fie fee? elie fenlesl ro 


=0 


Here (3) denotes a minor of the Jacobian oe of the transformation (2.1), 
B 


and therefore 
x 


OX = a a 
eel Oyk” = 0% 2* 


6; being the Kronecker delta, whence by differentiating with respect to e and 
then putting e =0 we have 


e ifs fay spe 5h + 0% Oy bx” = 0% 0, Ox" 


Thus, finally, by (2.4) we get (2.6). 
The functional 


d a 
(2.8) owael a f axt(x Bi} 
T e=0 
is called the local variation of the action W given by (1.20), and 
(2.9) AW = ae f dei (eB)} 


where 7 denotes the region obtained from t by (2.1), is called the total variation 
of the action. We have’? 


(2.10) AW = f dc { = 


sz OP" + 0,(L 82%}. 

So far, the variations of the field (x), both local and total, have been 
arbitrary infinitesimally small functions belonging to the whole space {. Now 
we impose on them the following conditions: 


10 Cf., e.g., COURANT-HILBERT [1930], vol. 1, Ch. IV, § 11, p. 8. 
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4° For every hypersurface o* for which dt, 6**=0, 


(2.411) Afar p= Os 
2° The variation 6% shall satisfy the equation 
(2312) 0, Op* =0. 


The variations 6)/% satisfying these conditions will be called the local hydro- 
mechanical variations of the field p. According to (2.4) we write 


(2.13) Ayb% = doh" + p6%- 5.x? 
and call A,p% the total hydromechanical variation of the field /. 


Theorem 1. All local hydromechanical variations are of the form 


(2.14) by b* = 05 (p* bx* — p* 6x"), 
where the dx* are arbitrary infinitesimal functions belonging to the space $. 


Proof. In order to prove that the condition (2.13) is necessary, let us sub- 
stitute formula (2.7) in equation (2.11): 


A f dtyp% = f dtq {Syp* — 8p (f° 6x* — p% 6x") + Og6"- bx} =0. 
From rhe conition 1° it follows that dt, 6x*=0 implies 
dt, {Oo b* _ Og (p° 6 x” — p* dx?)} =i) 
Oo b* — 0, (p" bx* — p* bx") = K 6 x%, 


thus we have 


kK being some factor of proportionality. But the condition 2° requires that 
0, (KK 6 x*) =0 for all d6x%*%, whence K =0, and (2.14) follows. The condition (2.14) 
is also sufficient. In fact, from (2.14) it follows that (2.12) holds identically, and 
from (2.7) that the condition 1° is also satisfied, which completes the proof. 

The formula (2.14) defines an infinitesimal group of transformations of the 
vector field #(x), depending on arbitrary functions 6 x* (x). 

We make a further comment on the conditions defining the hydromechanical 
variations. The condition 1° is an appropriate counterpart to the condition 
(1.13). In fact, in the particular case when CtcedV Ae 0y P12 ee 
condition 1° requires that when @ dV 6¢=0, and hence 6¢=0, (1.9) must hold. 
From (2.14) it follows also that then (1.10) holds. But this is only a special case 
of the general formulae (2.11) and (2.14), véz., that obtained by putting « =0, 
6x®°=6t=0. In the general case, the condition 1° requires that for variations 
6x* tangential to the hypersurface o* the total variation of the complete matter- 
flow contained on o* shall vanish, corresponding to (1.13). 

The condition 2°, which is equivalent to 6 (0,£*) =0, requires no local variation 
of the source density; this corresponds to (4.14), being a generalization of it 
appropriate to continuous media. It should be stressed that the condition 2° 
makes no use of the continuity equation 


(2.15) 0,~% = 0, 


which is thus superfluous in deducing equations of motion from the variational 
principle. 
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It may be noted also that the concept of hydromechanical variation can be 
introduced also if X, is a Riemannian space. It suffices to write the symbol of 
the covariant derivative instead of @,, as can be verified by calculation". 

2. After these preparations we formulate the variational principle: 

For all 6x* vanishing on the boundary of the region t, 


(2.16) AW = Af dz tia,p(x)) = 0, 


provided that the 6yp* are hydromechanical variations. 


We now deduce the equations of motion from this principle. By the formulae 
(2.10) and (2.14) we have identically in all variables 


AW= facto 03 (f° 0.x* — p*.x*) + 0, (L bat) 

(2.17) ; 

= f dx d,(TP dx") — f dry, 6x", 
where : 

Beets oL Bl y OL 
(2.18) ibs p apt On (L p a 
and 
rege OL\ 4 (ek 

(2.19) Ya P {2( 592) an (Saf 


Since the 6x* vanish on the boundary of the region t, the first integral on the 
right-hand side in (2.17) vanishes, and since the 6x* are arbitrary in the interior 
of t, we obtain from (2.16) the following scheme of equations of motion: 


(2.20) Wy = 0, 
which can be transformed into 

erie 1s OL \ ol 
ee p on (oa) | 2,(L p Pole ae 


where the sign 0L/0x?%, denotes the explicit partial derivative, 7.e., 0OL(y, p)/0* 
Loge 
The four equations of motion (2.20) are linearly dependent, since identically 


(222) PVa 0: 


(This fact will be discussed in Section IIT, 3.) Thus the four functions p* describing 
the flow satisfy only three independent equations. 

The same equations can be also obtained by using the local variation 6 instead 
of the total A, 7.e., by requiring that bW=0. 

In the case when the Lagrangean takes the usual form (1.21), we get from 
(2.20) for « =0 


= (-#?) +#grad (5?) +8 (5 grad P + grad U) == () 


11 A formula analogous to (2.14) is used by Fock [1955] pp. 214—218, particularly 
formulae (47.35) and (47.36), but deduced in a different way. 
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and: fon ==? —==4, 25.3 


ab) 


He + grad (5 8") — Tx rot T+ » grad P + grad U =0 


where 


P 


=o =-—é 
© To 
is the pressure. When « =0, the equations (2.24) give the same Bernoulli equa- 
tion as above, and for « =1, 2,3 the usual Euler equations: 
+ (v grad) v + ? grad P + gradU =0. 

All these equations of motion have been deduced under the assumption that 
the functions #*% are continuous in the whole region t. However, if there exist 
some hypersurfaces on which the #*% cease to be continuous, jump conditions 
for them can also be deduced from the variational principle. For simplicity, 
but without any essential loss of generality, let us assume that in t there exists 
such a single hypersurface o* on which the functions #* are not continuous but 
approach on each side finite and continuously differentiable limits. Let 


(2.23) F(x) =0 


be the equation of this hypersurface, and let its unit normal vector yg, which 
is proportional to the derivative 0,/, be continuous on o*. We assume that the 
hypersurface o* divides the region t into two parts t and t*, and that the normal 
vector is directed from the part t into the part t*. Denoting by dtz and dt; 
the elements of hypersurface directed inside t and t*, respectively, and by do‘ 
their common absolute value, we have 


dt; =—vgdo", dt; =v,do". 


Consider the expression 7? given by (2.18). Let Te denote the limit of 7% 


+ 
if we approach the hypersurface o* from within t, and let 7? be the corre- 
sponding limit from within t*. Then equation (2.16) becomes 


AfdtL+AfdtL =0, 


and from (2.17) it follows that inside both regions t and t* the equations of 
motion (2.20) hold. But on o* 


f do* v5 (T? — T?) 6x =0 


for all 6x*. Thus, denoting the jump of a function on the hypersurface o* by 
square brackets, e.g. 


ob 
[1] = 12 — Ty 
we obtain the following scheme of jump conditions: 


(2.24) isle hea O. 
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By substituting from (2.18) we get 


p OL 
Op* 


since ¢, F is proportional to »,. 


| pF IZ Fa 0,F =0, 


(2.25) p 


It is worth noticing that these jump conditions hold also when there is a 
distribution of sources on the hypersurface o*. In the special case when there 
are no such sources, then 


(2.26) [p°] a, F =0, 
and hence, writing 
fe g pb OL =, 
(2227) L—# ape +> 
we have from (2.26) 
rece oL 
(2.28) [P] 0, F = — (p? a F) ae: 


The expression P given by (2.27) is the generalized pressure, as in the case when 
the Lagrangean takes the form (1.21) P reduces to the usual pressure. The 
jump-conditions (2.28) then take the usual form also. In fact, let 7 be the space- 
like unit normal vector to the discontinuity surface, 7.¢., let 7 =grad F/| grad F, 
where “‘grad”’ denotes the spatial gradient. Then 


is the speed of displacement of the surface of discontinuity F(é, x1, x7, x8) = 0, 
and equations (2.28) become™ 


o(G—»,) |£9?— 4] =o, (PI, 
o(G —»,) [8] =#(P). 


III. Conservation laws 

1. In the calculus of variations the theorems of Emmy NoETHER™ describe 
a relationship between the invariance, with respect to given infinitesimal trans- 
formations, of the functional considered in a variational principle, and some 
identities satisfied by the usual Euler expressions of this functional. The essential 
idea in NOETHER’S theorems can be adapted in the variational principles of 
hydromechanics considered here, although the hydromechanical variation differs 
from the usual one. 

There are two kinds of Noether theorems. In one, the transformation is 
supposed to depend on scalar parameters; in the other on functions. Accordingly 
we shall formulate two theorems corresponding to the hydromechanical variation. 
In the application of the first theorem, the transformation of the space X, will 
be taken as the Galilean group. This will restrict the form of the Lagrangean 


12 Cf, LICHTENSTEIN [1929], p. 370. As mentioned by LICHTENSTEIN, ZEMPLEN 
[1905] was the first to obtained the jump conditions from a variational principle, 
using, however, Lagrangean variables. See also TauB [1949], p. 149. 

13 Cf, NOETHER [1918], p. 258—276. 
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and lead to conservation laws of matter, energy, impulse and angular momentum. 
In the application of the second theorem, the transformation of the space x4 will 
be taken as one for which the hydromechanical variation vanishes identically. 
This will lead to a generalization of the conservation laws for vorticity. Before 
stating and proving the corresponding theorems, we give some definitions. 

Let 
(3.1) W = fdr Ele h(a) 


be a functional defined in the space ‘8, and let t be any given region contained 
in X,. Consider the infinitesimal transformation of the space X, 


(3.2) R= a PAS, 


where the 4x* are functions of x, # and their derivatives ¢, p up to a given order. 
Consider, further, the infinitesimal transformations of the functions (x) given 
by the formula 


(3.3) P* = p* + 05 (f° Ax* — p* Ax?) + Ogp*- Ax? 


and corresponding to the transformations (3.2) of the x. By means of (3.2) and 
(3.3) the functional is transformed into 


(3.4) W=W-+A,W = fdrL(z,p(%)), 

where 7 is the transformed region t. By Siting 

(3.5) Sob" = Op (bP Ax*— pr dx’), 

we have from (3.3), (3.4) and (2.10) 

(3.6) AyW = [arf oe bo P* + 0,(L Ax”). 
If there exists a vector C* such that 

(3.7) A.W = f dx 8,C* 


identically in t, then we call the functional “hydromechanically invariant up 
to a divergence’’, or “‘div-invariant”’, with respect to the transformations (3.2) 
and (3.3)'. If, in particular, C*=0, so that 4,W=0, then W is said to be abso- 
lutely invariant with respect to (3.2) and (3.3). Finally, we call the expressions 
Yo given by (2.19) the hydromechanical Euler expressions. 


2. Theorem 2. If the functional W is div-invariant with respect to the trans- 
formations (3.2) depending essentially on k arbitrary parameters, then exactly k 
linearly independent linear forms of the Euler expressions Wo, ave divergences. 

The proof is based on the following identity 


(3.8) AW = J drag (Te Ax*) — yp, Ax}, 


“The notion of div-invariance was introduced by E. NortHER; cf. BESSEL- 
HAGEN [1921], p. 261. 
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given by (2.17) and (2.18). By the hypothesis of the theorem, the transforma- 
tions (3.2) are of the form 


(3.9) AE it ee 


where ¢” are infinitesimal scalar parameters, and /*%, some given functions. 
Besides, 


(3.10) Oe EC, 
where the C%, also are given functions. The div-invariance of W is identical in 


the parameters «”; therefore, substituting (3.10) into (3.7) and (3.9) into (3.8), 
we have from (3.7) 


e” {dx O,Ch, = & f dr {O,(TP1%) — yp, Ie}, 
identically in «” and in t. Hence 
(3.11) en eae EA fo oem Oe ORD ee? 


which completes the proof. 


During the motion we have 


(3.12) Ya = 0. 
It follows from (3.11) that the functions * satisfy the equations 
(3.13) Onl T? Uy, = Ci) = 0 


which are called in physics “‘conservation laws’’, usually being written in integral 
form. Namely, let o* be any closed hypersurface consisting of two parts o} 
and o;,. Integrating (3.13) over the region bounded by o* and transforming the 
integral by Gauss’ theorem, we obtain an integral form of the conservation laws: 
(3.14) J deg (TP ls, — Ch) = J dt (Te by, — Ch)» 


am 


. 


oF Orr 

which asserts that during the motion (3.12) these integrals are conserved. From 
the general scheme of conservation laws given by (3.13), particular forms can 
be deduced in accordance with the particular kind of transformation (3.9) assumed. 
We shall confine ourselves in applying formula (3.13) to the case when (3.9) is 
the group of Galilean transformation, basic for classical hydromechanics, although 
the same method can be adapted to other transformations such as the Lorentzian 
group. Suppose, then, that the action is absolutely invariant with respect to 
the Galilean transformation 


(3.15) Ae ie, 0 mn cc, B= 0,152, 3, 
in which the infinitesimal parameters a* and aj} satisfy the conditions 
(3.16) ag—a=0, a@+a=0, 1,7 =1,2,3. 


In view of the absolute invariance assumed, in (3.7) we have C%=0, and, therefore, 
in (3.6) we have 4)W=0, identically in the parameters a and in the regions tT. 
Putting (3.15) and (3.16) into (3.5) to (3.7), we get 
ts eee aly sent ign 2 

“Opi ie api ap* | OF 


One, 


(3.17) p 0, %p° 
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From the first of these equations it follows that the Lagrangean L can depend 
only on the sum (p!)?-+ (p2)?+ (p3)? and on f°. Thus the invariance of the action 
with respect to Galilean transformations restricts the form of the Lagrangean). 


From the second equation (3.17) the continuity equation d, ~p®=0 follows. In 
fact, d,p°+0 would imply, by (3.17), ae oe = '/p', which is impossible, 
as the expression on the left-hand side depends on the sum (f1)?+ (p?)?-+ (p%)?, 
while the right-hand side is independent of it. 


Other conservation laws are given by (3.13), which yields in this case 
(3.18) 6,7 =0, 817 =0, a%ME=—0 


with Mi = Tf Few T’x;. Substituting for T? the expression given by (2.18), we 
get the general scheme of conservation laws for barotropic flows of an inviscid 
fluid, viz., that of energy, impulse, and angular momentum, respectively. If L 
takes the particular form (1.21), these laws become the usual ones. 


3. Theorem 3. I} the functional W is div-invariant with respect to the trans- 
formations 


(3.19) Ax* = AS (p'] = AS gy + ASM Oy, Go + ASH Oy aD 
depending essentially on rv arbitrary functions g*, s=1,2,...,7, and their deriva- 
tives up to a given order q, the coefficients A being given functions, then there exist 


exactly v linearly independent identities between the Euler expressions wy, and their 
derivatives. 


Proof. By the hypothesis, in (3.7) we have 
(3.20) (CO Re [e*] = Ce yg a Ca" 0a, yp Ses Ce Oa ade Q’, 


where the C¥~ are given functions. By substituting (3.19) into (3.6) and (3.8) 
we have 


(3.21) fdr d,(I*[o']) = f dr{a (TZ AS (y'] — y, At [¢*}} 


identically in the gy and in the regions t. Let Ae be the operator adjoint to the 


operator A%, 7.e., 


(3.22) AS [p"] = Ab gp’ — On, (AR) + + (=) gag (ASQ). 
Integration of (3.21) by parts yields 


eo Jaca; [yal p= $ day (IP Le’) + TP AS lg) — 
eas $ dt, (AS pe gp + Atht Ya Oy, GP ee a Ata Aggy Oe P). 


As the functions g* are arbitrary, we select g* which vanish, along with their 
derivatives up to order g— 1, on the boundary o* of t. Therefore, from the identity 


** Apparently, E. & F. Cosserat [1909] were the first to consider the dependence 
of the form of the Lagrangean on the invariance properties with respect to a given 
group of infinitesimal transformations. Cf. op. cit., pp. 8—10, 68—71, 125—127. 
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(3.23) it follows that 
(3.24) A®Ty,] =0, 


and these are the identities for the Euler expressions y,. 


In order to apply Theorem 3 to the hydromechanical case, we have to know 
the explicit form of the transformation (3.19), which in general depends on the 
form of the Lagrangean assumed. However, one group of transformations (3.19) 
can be produced for which the action W is always div-invariant, whatever the 
form of the Lagrangean may be. This group consists of those 4x* for which 
09 P*=0 identically, 7.e., the hydromechanical variation vanishes, as follows 
immediately from (3.6) and (3.7) by putting C’=LAx*. Again, one particular 
case of this group is evident, viz., Ax*=p*q@ for arbitrary functions gy; thus 
y=1, A*=#*, g=1 in (3.19). In this case, Theorem 3, viz., formula (3.24), 
yields the identity p*y,—=0, showing that the four equations of motion deduced 
from the variational principle are linearly dependent, as already remarked in 
Section II. There we saw also that for the usual form of the Lagrangean (1.21) 
one of the four equations of motion is the energy equation. Now we see that this 
fact does not depend on the particular from of the Lagrangean. 


But the group 4x*=%q is only a special case of the more general group 
given by 

Theorem 4. The hydromechanical variation dyp*= 0,3 (p° Ax*— p* Ax?) vanishes 
if and only if 


(3.25) Axt = pp + “8 exh, 


pu, 


where p is an arbitrary scalar function, uz 1s an arbitrary vector, and p, is any 
vector satisfying the equations 


(3.26) (8, P, — 8491) = 0. 


Proof. First we prove that the condition is sufficient, 7.e., that for the 4x 
given by the theorem we have 6)f/*=0. Let us write 


(3.27) [oP = — fox = eth H 8, 
thus 
(3.28) 0, Pu _ Ci, P; == Doesnt e and Og jh 0e 


From (3.28) it follows that 
p’ Oris =0, 
which can be written as 


(3.29) pr fP7 He p fe ae p” ee Ti 


Substituting (3.25) into the expression for 69% yields, in virtue of (3.27), (3.29) 
and (3.28), 


bo b* = Op (9° Ax — p* Aa®) = 0p {— (pp + 2 /?”)} 


== PY Uy jap\ __ ii 
an (Fatt) af? =0, qed. 
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Now we prove that the condition given by the theorem is necessary, 1.0, 
that (3.25), subject to the restriction (3.26), gives all solutions of the equation 
69i*=0. Let us write 


(3.30) PRAT Reel eae 


Thus /** satisfies the equations 
Oe, ie =0 
and 
BPP pe ho eee 
The first of these equations implies that 


(3.31) eae a, Pu» 


where , is a vector, and the second one yields (3.26). Let u, be any vector. 
From (3.30) we obtain 
— p*u, Ax’ + pu, At = f°? ug, 
and hence, writing 
Uy AX? _ 


nie tna 


we get 


eau a, Pu 


ee oes Up eae piled Up 
Ax*=>p CAMs! =f Ol ree 


what, together with the equations (3.26) already proved, completes the proof. 
For later use we remark that the expression p°Ax*— p*A.x* does not depend 
on the particular choice of the arbitrary vector #;. This follows from the for- 


mulae (3.30) and (3.31). Therefore, all expressions 4x* calculated for various ug 
by (3.25) differ by a term proportional to the vector £* only. 


4. Now we deduce from Theorems} and 4 a generalization of the vorticity 
theorems of hydromechanics as conservation laws. It would seem to be simplest 
to express the 4x* satisfying Theorem 4, 7.e. as given by (3.25), by means of 
independent arbitrary functions without any side condition like (3.26), and put 
them into (3.8). However, such a method is rather cumbersome. Therefore, 
we shall adapt the idea of Theorem 3 to this case, when instead of (3.19) we 
have (3.25) with the side condition (3.26). 


Thus by (3.6) and (3.8) we have for 6) #*=0 the identity 
J dt py Ax = f dv n{(T8 — OBL) Ax. 


We transform the integral on the right-hand side by Gauss’ formula into a hyper- 
surface-integral over the boundary o* of the region t considered. On o* we take 
Y =,=0 and the vector uz as normal to o*. Then 4x*=0 on o*, because in the 
expression (3.25) the derivatives are taken in directions tangential to the hyper- 
surface o* only. E.g., in the expression for 4x! the derivatives of the component 
Po are taken in the direction of the vector w, e'”?° perpendicular to u,, which has 
been assumed perpendicular to o*. So, by using (2.22), we get 


(3.32) [ards = [dry te oa, 9, =0, 
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provided that the side conditions (3.26) hold. These side conditions we take 
into account by means of Lagrangean multipliers 2”. Thus (3.32) yields 


uU o4 L y;} A L ‘. —— 
far (vezite pant O p" —Q! ?”) a, Py, = 0, 


where, now, the functions g, are arbitrary, provided that they vanish on the 
boundary o* of t. Integrating by parts in the last identity leads to the following 
conclusion: There exist vectors 2” such that 


(3.33) 0, re ex Any.) = 0,(Q" p* — Q’ p"). 


These are the identities corresponding to (3.24) in Theorem 3. 

Now we shall show that these equations are really generalizations of the 
vorticity-conservation laws in hydromechanics. To this end we remark first that 
the functions £2“ do not depend on the particular choice of the vectors uz, because 
the integral on the left-hand side of (3.32) does not depend on wg, as follows from 
the final remark in §3 and from (2.22). Thus in (3.33) the vector u, can be 
taken as completely arbitrary. Secondly, we remark that one of the Lagrangean 
multiplier Q“ can also be chosen arbitrarily, since the side conditions (3.26) are 
linearly dependent in view of the identity 


bp’ p" (8, Pu a On Pi) =. 
$= 0) 172930) The identity (3.33) gives for w0 
(3.34) 0;(Q'p°) = 0; 
for “=1 
(3-35) ($5) — 5 (-$t] = G02" p%) — 3, (Q%p! — O1p%) + a (Q'p* — 29); 


and analogously for u = 2, 3. 


Therefore, put %) =1, 


The meaning of these identities will be clearer if we introduce the usual 
three-dimensional symbolism, denoting the spatial vectors (p°Q!, f°Q?, £°Q3), 


(Wr, Ye, Ys) and Ga de r) by Q, w and v, respectively. Then (3.34) takes the 


p?’ Pp 
form 
(3.36) diy 2 ="6; 
and (3.35) becomes 
p90 i tiapls 

(3.37) rot 5 ee tot (Qx iv), 
or, by use of a known formula, 

OO ices ae Rn 2 
(3.38) rot a ae + (% grad) Q — (Q grad) ¥+Q divi. 
During the motion we have y =0, and (3.38) reduces to ZORAWSKI'S criterion!6, 
(3.39) as + (@ grad) Q — (Qgrad)# + Qdivé =0, 


for conservation of the vector-lines of the field 2 and the strength of its vector- 
tubes. 


16 ZORAWSKI [1900]. Cf. also TRUESDELL [1954], pp. 55—57. 
Arch. Rational Mech. Anal., Vol. 2 28 
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It remains only to show that in (3.38) and (3.39) the vector 2 can be taken 
as a curl of the vector OL|Op with the components (dL/0p1, 6L/dp?, OL/0p%), 
Lil Q =rot dL|0p. In fact, calculating the vector ~/p® occurring in (3.38) by 
means of (2.19) we get 

OL 


y _ o (OL gee ex YOU 
5 = 5 (25) eshte th 


IE, oan 
and by using this formula we show that the vector rot — satisfies (3.38) and 


(3.39). If the Lagrangean takes the usual form (1.21), then OL/ap= v, so that Q 


becomes the usual curl of the velocity 2 =rot v and the identity (3.38) gives the 


usual Helmholtz equation. 
However, formula (3.39) shows that in addition to rot 6L/ép there exist 


infinitely many other vectors 2 whose vector-lines and intensities are conserved 
during the motion y=01; indeed, it would be sufficient that rot py =0. 
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AA Method of Solution 
of the Equations of Magnetohydrodynamics 


JUDITH BLANKFIELD & G. C. MCVITTIE 


1. Introduction 


The equations of classical magnetohydrodynamics for an infinitely conducting 
medium fall into two groups. The first group contains the three equations of 
motion and the equation of continuity, the second group contains equations which 
involve only the magnetic field and the velocity of the gas. The first group of 
equations can be thrown into the form of the vanishing of the vectorial divergence 
of an energy-tensor (see Eqns. (2.14), (2.15) and (2.19) below). This group can 
be solved in terms of indeterminate functions in a variety of ways, for example, 
the very general type of solution of equations of this kind found by TRUESDELL 
(1957) in his study of a problem in elasticity theory could be adapted to the 
present purpose. Alternatively, the theory already developed by us (BLANK- 
FIELD & McVITTIE 1959, hereafter referred to as BM), which depends on de- 
generating Einstein’s equations of general relativity, may be employed. In this 
method the indeterminate functions are known to be degenerate forms of the 
generalized potentials of relativity theory, and one of them, indeed, is the ordinary 
gravitational potential of Newtonian theory (see the first of Eqn. (2.21) below). 
But readers who are less interested in the origins of the solutions than in their 
application to specific physical situations may pass directly to equations (2.21) 
and verify by direct substitution that the expressions there given do, in fact, 
satisfy the equation of continuity (2.11) and the three equations of motion (2.14). 


In what follows the motion of the gas is regarded as controlled by its pressure- 
gradient and the magnetic field. The specific illustration which is worked out 
in detail refers to one-dimensional magnetohydrodynamics in which the de- 
pendent variables involve the time and one linear coordinate only. The physical 
situation is that of a slab of gas, of infinite lateral extent, with parallel plane 
boundaries. It is shown that in adiabatic motion the slab can expand until it 
obtains a finite width and an equilibrium configuration. 

The same notation is used as in BM: Greek indices range over the values 
4,1, 2,3, the index 4 being associated with the time, the others with the space 
coordinates. Latin indices take on the values 1, 2,3 only and /mn always 
denotes a cyclic permutation of 123. A comma denotes partial differentiation 
with respect to the coordinate(s) indicated by the subscript(s) following the 
comma. 

28* 
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2. The Equations of Magnetohydrodynamics 
The equations describing the flow of an infinitely conducting non-viscous 
fluid are (COWLING 1957) 


oH 20. (2.01) 
V,H+Vx{HxUl|=0, (2.02) 
oV,U +0(U-V)U+Vp+uHx[V x HW) =0, (2.03) 
Oca Ve 0) =i) (2.04) 
where 
J=VXH, (2.05) 
E=y(Hx0U), (2.06) 


and U is the gas-velocity, 9 the density, # the pressure, H the magnetic field 
vector, E the electric field vector, J the current per unit area, w the permeability, 
and V, the partial derivative with respect to the time, mks units being employed. 
For comparison with BM it is necessary to convert these equations to a dimension- 
less form. Let J, t, V be constants with the physical dimensions of length, time, 
and gravitational potential (mass/length), respectively, in mks units, and let h 
be a constant such that 


Ue me 


— (2.07) 


Let the vector H have the three components # H; where the H; are dimensionless, 
and let U have the components U,//t where the U; are dimensionless. Further, 
a dimensionless density can be obtained from @ by multiplying it by /?/V, anda 
dimensionless pressure from p by multiplying it by #/V. Using 0, p to denote 
these dimensionless variables, (2.01) to (2.04) may now be written 


3 
2 H;,;=0, (2.08) 


A, 4 me (H, U, m ree Uy) as vr (i, U, = H; oe = 0, (2.09) 


Q Un, i 03, U,; ste Pi =e re (Ain 1 —- fie) =F 7, (A, , ae H,,1) aa O, (2.10) 
j= 
3 
@at Lilo U0. (2.41) 
j= 


If equation (2.11) is multiplied by U; and added to equations (2.10), the following 
equations are equivalent to (2.10) 


(e U), +3 (e UU; Vp Dai i dd etna Hy == TT El on ett is ee Oe (212) 


By engine and subtracting the expression H, H,, and adding the expression 


— H, pa ; which, by (2. 08), is zero, (2.12) may be written as follows 
i] 
3 
(0U),4+d (oU,U;+ 6; ),; Bi 2h Hi — HH, ;— HH; ) =0, (2.43) 
j= 


j=1 1 
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or 


a 2 2 2 
(eU)s+> Je U,U; — HH, + 61; (6 4 abet at =0, (2.14) 
j=l 


7] 
which replace the equations of motion (2.10). 


A formal solution of the equations (2.11) and (2.14) can be arrived at by 
analogy from the case when the magnetic field is absent. Under this condition, 
(2.11) is unaltered but (2.14) is 


3 
(UN a leu, U; + 6;;)),; =0. (2215) 
i= 


The ten functions oe, 9 U;, @U,U;+6,; p, are the components of the Newtonian 
energy tensor, 7“”, of the fluid and the equations (2.11) and (2.15) are equivalent 
to the statement that the vectorial divergence of this tensor is zero. There is 
a corresponding formulation in general relativity of this problem, namely, if 
7" is the energy tensor of the fluid then, as is shown in BM, this tensor is 
related to the metrical and Ricci tensors (in dimensionless variables) by Einstein’s 
equations 


v AV 1 Lv 
eel | l= RACER (2.16) 
where 
8uaGV 
eh aaa, (27) 
cee. (2.18) 


Cc 


G is the constant of gravitation and ¢ the velocity of light. It is also known 
that the vectorial divergence of the right sides of (2.16) vanishes. Therefore, if 
we calculate the right sides of (2.16) in the Newtonian approximation we obtain 
expressions for the components of 7”“” in terms of the functions involved in the 
Newtonian forms of k“’— }g""R. But now the magnetohydrodynamic equations 
(2.11) and (2.14) differ from the hydrodynamic ones (2.11) and (2.15) only in 
that the components of the Newtonian energy tensor are 


Tt — 0, 

T*'= 9U;, 

Tie 6.0 a Hee 

T= oe Upe? HSE, + A). 


(2.19) 


Thus for the corresponding general relativity problem, a new energy-tensor 7,” 
would be needed in (2.16) whereas the right sides could be retained unaltered. 
Degenerating to the Newtonian case, the components of the energy-tensor would 
be given by (2.19) and these would be equated to the degenerate forms of 
R“’—i¢"*R. In BM, the right hand sides of Einstein’s equations (2.16) were 
worked out for the metric 


do? = [1 —n(p + 22° y,)| (464)? — 


3 3 . 
AVC! MN aaah CURE ae 
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where o is the dimensionless arc length, &“ dimensionless coordinates, and y, 
y,, and ¢; are eight dimensionless indeterminate functions of the &”. With the 
energy tensor given by (2.19), the Newtonian approximation to Einstein’s equa- 
tions (2.16) is: 


Cr iy, 
3 
0 U, =yath(V%C oe dba , 
Y U; UZ a 7, A, i (Wa sr Dahtan - z ay Gita) , (2.21) 


e UP +P +3 (— AP + Hm + Hn) 
= Pian Sade Cay (Wn Wa) ips en Va) a 


These can be regarded as solutions of the equations (2.11) and (2.14) for an 
arbitrary choice of the eight functions y, Wyo;, as can be verified by direct 
substitution. However these eight functions must, of course, be determined so 
that the magnetic field and the velocity satisfy (2.08) and (2.09). It may also 
turn out that the form of the energy-tensor adopted in (2.21) imposes inter- 
relations between the eight functions, called consistency relations (MCVITTIE 
1956). 

In one-dimensional motion, all dependent functions involve, by definition, 
one linear coordinate and the time. Let X be the dimensionless Newtonian linear 
coordinate, here designated by the index 1, and let T be the dimensionless New- 
tonian time, designated by the index 4. The equations (2.21) reduce to 


Oars Vans (2.22) 
oU, =a 
@U,=$ loi, (2.23) 
0 U3; =3 lair, 
oU, US hie 2,14> 
0 U,U;— H, H, =0, (2.24) 


@ U; U, — H, A, = — = Cs 145 


eUl+p+3(—Ait+ H}+ 2) =—yay, 
0 UZ p races AY H3) = = Pagttial + Ys— Va a: (2.25) 
@U3 p+ 3(Ai+ A H3) Sere aaa Chat t (po W4),11- 


Maxwell’s equations (2.08) and (2.09) are, in this one-dimensional case 


Ay y — 0, 
fH, ,4=0 
Tae (ear Ui) 4 =0 


As 4 — (Use, Uy) 4 =0 


(2.26) 


d 


’ 


and therefore the component, H,, of H in the direction of X must be a constant. 
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3. A Particular One-dimensional Solution 


The special case of the solution given by equations (2.22) to (2.26) in which 
the ¢; are taken to be zero will now be examined in more detail. The equations 
(2.22) to (2.25) are now 


Ce Ts (3.01) 
oU,=Yar 
OU = 05 (3.02) 
OCs): 
oU, U, = Hal, = 0, 
0 U,U, — H,H,=0, (3.03) 


o U; U,— H, H, =0, 
OU, Pig — a a) Yas 
0 Uy + p 3 (Hi H} H3) Yaa + (Ps Ya),11> (3.04) 
o U3 ot 3 (Hi = H; =F H3) = — Paat (Po Wa),11- 
From equations (3.02) and (3.03), it is easily seen that 


I 


U, + 0, U,=0, U; = 0, 
H,H, = H,H, = H,H,=0, 


so that the motion is parallel to the direction of X and at most one of the magnetic 
field components is non-zero. The non-trivial case is that in which the magnetic 
field is not constant, that is, when H,, the component perpendicular to the 
velocity, is the non-zero component. In this case, H, can be found from the 
third of equations (2.26) after a substitution has been made for U, in terms of w 
from the equations (3.01) and (3.02). Thus 


Wray a 
Hy«+(— Yur H,) =0, G0) 
or 


ieee 6 = eee 
Ue ae Hse ola (3.06) 


This equation can be solved* with the aid of the equation of continuity (2.11). 
Dividing equation (3.06) by H,, and the equation of continuity by @, gives 
Ay, ves aten Ay y Bt (sat | 
H, Yat A, Yair Be 
BL Ae VU OL (#1) : 
Q Vibe. 12 Pelt 
The difference between these is 
Ay 4 Q,4 ee 22} 
as ; ; : —| =O, OW 
A, Q Wii (3.07) 
which has the solution 


r. H,=0 Dy 3) Se Ply 1) , (3.08) 
* This solution was suggested by A. H. Taus. 
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where @ is an arbitrary function of y,. Regarding this as the solution for H, 
in terms of yp, the first of equations (3.04) may then be used to determine p. 
The last two of equations (3.04) merely serve to determine (y,—¥y,4)11 and 
(w3—a),11 When everything else is known. Thus, these last two equations are 
consistency relations and play no part in finding a solution of equations (3.01) 
to (3.04) 

To show how these steps can be carried out, a particular solution will be 
obtained for the case of linear waves, by which is meant that the velocity is 
linear in X, say 


ee 09 
al See Ya4 F ? (3 ) 
where F and G are arbitrary functions of 7. A prime will denote a derivative 
of a function with respect to its variable, in the case of F and G therefore, with 
respectito [ket 


x . 
Z=5+6; (3.10) 
then 
yi i= F(Z), (3.11) 
yp =FF(Z) +H, (3.12) 


where H is an arbitrary function of T, and ¥ is an arbitrary function of Z. 
From (3.08) 
H,=F(Z)/E, (3.13) 


where ¥ is an arbitrary function of Z. The pressure may now be found from the 
first of equations (3.04) as follows 
p=—F"F 4+ (ZF — (FG)") F'— HH" -1 FF? 
a 
=F"(Z F'— F)\ —(FG)" F'— A" —2 Z2IF?. Oy 
It is now necessary to introduce boundary conditions. Suppose that the 
material occupies an infinite slab bounded by planes perpendicular to the direction 


of X, and that the pressure and density are taken to be zero at the plane boundaries 
and also outside the slab. Since 


1 vr 
let the function ¥’’ vanish for two values of Z, say Z==0, and Z =¢,, so that 
these points determine the boundaries of the gas cloud. Assume 
FCO! 
Eid VANE A Oe 6) VA ae 


Define constants A, B;, and C,, for.z=41, 2 as follows 


F'(x;) = B,, 
F (a;) = C,, (3.47) 
I (o&;) =A 
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The condition that the pressure is to vanish at the boundaries where the density 
is also zero may be expressed by the following two equations 


0.= F(a, B;— C) —(FG)B;— A" pA 


a 2 F2? (3.18) 


the difference between which is 


PG, B, — a, By, — C, 4 C2) a ea) G)"" (By = B,), 
or (3.19) 
[F(G — 2) =O); 


where 
Ba mF = Bef (3.20) 
which has the solution 
Gee AB) (3.21) 


where « and f are arbirary constants of integration. The two equations (3.18) 
now yield 
‘eel 


Jah! eI 6, 5 ET? (3.22) 
where 
Gre By B.—% ae C,—B, Cy ; (3.23) 
Thus the expression for # is 
’ x! wre Albee 
p =F" {(Z — B)F'— F —C} + =. (3.24) 


In the external region, the magnetic field vector must satisfy Maxwell’s 
equations for empty space. The latter reduce, in the one-dimensional case, to 
a single equation for H,, namely 


Ay = 8 Ay aa, (3.25) 
where ¢=(I/t)/c as in Section 2. The general solution of equation (3.25) is 
H,=/(T + eX) +e(r Sex), (3-26) 


where / and g are arbitrary functions. The boundaries occur at X =(%,—G)F 
and X=(x,—G)F, and H, must be continuous there. Hence, using (3.13) and 
(3.17), the boundary conditions are 


PE + € (oy G) F) + g(T € (Oy G) F) = 


f(T + €(«,—G)F) +¢(T — e(,—G)F) = 


The general solutions of these equations for / and g in terms of F will not here 
be attempted. It is however possible to obtain a solution fairly simply in the 
symmetric case, that is, the case where the center of the gas cloud remains at 
rest at the origin. Since the boundaries are at (%,—G)f and (%,—G)F, the 
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center is at (2% —G) F. From (3.21), the center 1S. aye X= ("5% — B) F— 


« IT —f. The symmetric case is defined by 


n= Pp =O, (3.28) 
Bae (3.29) 
Banaeos (3.28) place a restriction on G, and equation (3.29) places a restriction 
on ¥. The restriction on FY is not an ee one to satisfy, as one could take 
F =cos Z where B =0 and «,= —a,=32. In the symmetric case the boundary 


conditions (3.27) are 
farang) ele eae, 
F 


‘(7 pI F) 4 g(T e= = F\= 


(3-30) 


The two conditions (3.30) become one condition if g is taken to be the same func- 
tion as f, thus 


Ln CAS Ci ae A 
Mie : F)+i(T ae F\=5, (3.34) 
and a series for the function / can be found as follows. Let 
f= De" ty, (3.32) 
N=0 
and expand each of the functions fy in a Taylor series 


fy(T£e%5% PF) = S(eeBS*R(T)) 1 I(T), 


n=0 


where /@ (7) means the m-th derivative of fy with respect to T. Hence 


f ai(T+e%S"F\+ (7 — 285A) 


5 OO 2wtn) (Ma laa i 
eee pa rn (aS 1 F(T) Gall (2m) (T) , 


and the functions /y may be found one at a time by equating powers of ¢ on both 
sides of equation (3.33). For °, 


(3.33) 


4 
eae 
for ¢?, ; 
h+5 (35% FY fi =), 
A feo 2 
memes Us. 
and so on. Thus outside the gas cloud, 
Sf (Le) eee eX), (3.34) 
where A A \2/ 4 \" 
ha 2x Oy — o 4 
~ feed gay yy, 039 


for the symmetric case. 
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To proceed further with the solution within the gas cloud, one can impose 
the additional condition that the rate of change of entropy, 


dS d 
GHE as Gk [In p yin a , (3.36) 


be zero so that the motion is adiabatic. The calculation of this quantity is 
simplified in the case of linear waves because the derivative following the motion 
of any function of Z alone is zero. Thus 


ie ee 
Figen F° 
PF (3.37) 
5 FF" (Z—B)F'—F—C} (42?) 
3 
Dit Safe i ; 
FF’ {(Z—B) FF —C} + a (A2— IF?) 
so that 
| wt eS vt gel K ie Y ae 
we OE ty F(Z B) F¥'—F—C} ac 4) (a2 J?) 
FAD eae ; 39 141338) 


F'{(Z—B) F' —F—C} + ry (AMP? 


The adiabatic condition, then, is 


wtr 1g aA y Fe by 
Fy F(Z — BF —F —C}=F (tL) (2-99), 6.39) 
or, assuming that F is not constant, 
B ety F*| (1 = ¥) (A2— 2) 
F’ —F (2B) Fa Be aGy G40) 


The left hand side of (3.40) is a function of T alone, while the right hand side 
is a function of Z alone. Therefore, both must be constant, say C*. When the 
right hand side of (3.40) is set equal to C*, assuming y= 2, one may solve for #? 
to obtain 

C* 


wi 
2 


$2 = AP Za) = P= GC) (3.41) 


1— 


When the left hand side of (3.40) is set equal to C*, the following equation is 


obtained 
(GEA JER 


gc yf’ J RY — Fe i (3.42) 
which can be integrated to give 
— , C* 


if y+1. Equation (3.43) can be integrated if y= 2, after multiplication by the 
integrating factor, 2F’~* F’, to give 


BP (Pi ee ZOLA, (3.44) 
V2 Vial 
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The expression for .£? found in equation (3.41) and the above equations for f 
may be used to obtain the more compact expression for the pressure 


p =A, Yo" (2 —B)F'—-F —C}. (3.45) 


Since the pressure is to be positive in the interval a;<Z<«,, it must be 
shown that the constant A, is positive, which can be done as follows. Clearly, 
the quantity }(y—1)F” is positive since F is positive and 1<y<2. Let 

(2) =C-—bF —A— C. 
The constants B and C are such that /(Z) vanishes at the boundaries of the gas 
cloud. Differentiating with respect to Z, 
['(Z) =(4— B) F(Z). 
It has already been assumed in (3.16) that 


F(Z) 20 = Si ae, 
so that 
((Byi20F “aa Se, 


$(Z) Oe 3B aaa 


Thus /(Z) is an increasing function of Z for 4,<Z<B, and a decreasing one for 
B<Z<a,, or {(Z) is positive in the interval «,<Z<«,. Hence, in order to 
have a non-negative pressure, it must be the case that 


Aa 0: (3.46) 


From the equations of motion (2.10), the force per unit mass including the 
pressure gradient is seen to be 


1 
F; ey rt [e.+ fa Pie bi vv Fs a) a TAS Ey an H,, 1) | , (3.47) 


and therefore in our case, by 


ee Vee (4 Ved ey Gy pe 
1 = ( ){ As 5 F ee 
Rane 


’ 


(3.48) 


If the expression A, ? FY (C*F-*)/(2—y) is positive, then the force, F, 
is everywhere directed outward, #.e., 

T= OF 5” oe ee 

SOS (BZ aan 
If, on the other hand, the expression A, — FE ?— (C*F*)/(2—y) is negative, 
then the force is everywhere directed inward, 7.e., 

Lis 05" is ee Be 

Eizs0s> (Bf 
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Thus, if the constant C* is negative, the force must be directed outward and the 
gas cloud must expand. But if 


then it is possible for the gas cloud to hold together. 


For arbitrary y the equation (3.44) is not soluble in terms of elementary func- 
tions. However if y = 3 it can be solved, and this case will now be examined in 
greater detail. Equation (3.44) is now 

Macs he iAd 


(F')2 ee + 4A,. (3.50) 


In order that F be a real function of 7, the right hand side of (3.50) must be 
non-negative, which will be the case only when 


Frte A, ( A, at 


8C* 8C* B, 
or 
Bie de AL etal 
B'> sce 4 (Gea eh 
Thus F will always be real if 
eee ee 
(scx) — Gt So, 
or 
A 2 
(Gh pe 
eye (scx) 
In particular, consider the case where 
Ave 
Equation (3.50) becomes, in this case, 
poli PSOE ee: ieee Cae 4 4 
é F rata oo ace) | ee F 8C* |’ 052) 
which has the solution 
jy ee te RE Yipee Wb Coie . SC al 
da A, ronal ee In (1 ace FH noe (>F nore i eo 
As T tends to + oo, F tends to the finite value ay, from below. Since the 


width of the gas cloud is directly proportional to F, the cloud expands, but 
only to a finite width. Thus a solution of the equations of magnetohydrodynamics 
has been obtained describing the motion of a gas cloud with infinite plane faces 
normal to the direction of motion, the forces being the pressure gradient and a 
magnetic force parallel to the faces. The pressure and density vanish on the 
plane faces, the motion is adiabatic and the ratio of specific heats has the value $. 
The gas cloud expands, but not beyond a certain finite width. Ultimately, the 
binding effect of the magnetic field balances the disruptive effect of the pressure 
gradient, and the cloud approaches an equilibrium configuration. 


This work is based upon a dissertation submitted by JupITH BLANKFIELD in 
partial fulfillment of the requirements for the Ph. D. degree at the University of 
Illinois. 
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On the Convergence 
of the Rayleigh Quotient Iteration for the Computation 
of the Characteristic Roots and Vectors. II 


A. M. OSTROWSKI 


In what follows I add some observations to Part I of this paper*, which the 
reader is presumed to have at hand. The enumeration of sections and formulae 
continues that of Part I. 


18. In the paper by CRANDALL referred to in §17 the cubic character of the 
convergence is proved in the following sense for the sequence of vectors &, obtained 
by the rule of § 9: 

If &,—>n, and if for a sequence of numbers ¢, with €,—>0, 


(58) E,— 9 =O (&) (k +0), 
then 
(59) E41 — 7 =0(8}) (k +0). 


However, from this result it does not follow that a, ,,;—o¢=O((A,—o)%), and 
still less that the asymptotic formula (46) holds. 

On the other hand, the exact relation (46) presents not only theoretical but 
also practical interest because by use of this relation the convergence can be 
hastened still more without much additional work. 


19. From (46) follows 


An — Ar +1 = Ay, — Oo — (Ap44—0) SEC) i 2 
A, — 6 A, —o : r (Ae 0) ); 


(60) 


Me Oey wien Ageia OTA, Aga as 
Therefore, by applying (46) again and replacing k by k—1 we see that 


An — Agar >y 
(Ap—1— Api)? 


, 


This paper was prepared in part under a National Bureau of Standards contract 
with the American University, Washington, D.C. I am indebted for discussions to 
Dr. E. V. HAYNSWORTH. 

* Arch. Rational Mech. Anal. 1, 233—241 (1958). 
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and by comparing this with (46) and again using (60), since y=-0, 


Avi 9 Aa Ait 
(Ap—Ap+1)® (Ap—1— Agtr)? ’ 
(An—Ants)* 


SRL ae aeae eee 


We see that replacing A,,, by 
' (An+1—An)* 
(61) Ana aa Ansa eS ae , 
gives an even better approximation to o than /,,,. Using this remark system- 
atically, we shall apply (61) at every second step. Starting with £,=7 we compute 


ho, ey; Ay, &s, A» 


by (28) and (29) and replace A, by Aj, using (61). Then we start with A, and & 
and proceed in the same manner. The general double step starts with Aj,» 
and &,,,-9. Lhen we compute é,,,_1, Aamir €em> 4am USing (28) and (29) and 
then again replace A,,, by A;,, using (61). As a result of a discussion published 
elsewhere*, it is expected that the gain in precision due to the above improve- 
ment is about 114%. 


20. In computing practice the cubic character of the convergence appears 
only at the second step. At the first step the error is approximately only squared. 
Indeed, we have from (13) for k =0 


; ol =o 
ot ee ae 
e. [— ’ ae 3 ¥ 
ae) p+) ease Pu 


But from this it follows, if ~2j>>o, that 


(62) A,— 0 as 1 yy: Pu 


(Ay—o)? Pp phe 


In order to discuss in an analogous manner the order of A,—o, we rewrite 
(34) for k =1, making use of (33): 


A. 0 == 
o a) p | m . a , 
(o— Ap)? (o—A,)? TER Mi. 1 
or 
p bip =O é 
x Ay— eee 2 jot nt ; in 
(Ay—a)? (A, —o)?® “ Pu 


p+(o—A)*(o= A)? SY 


re (6,— Ay)? (6, a A)? 


Now if we let 4) approach o, we see that J,—o, A,—>o, and it follows that 


sth Oe eg ve ae, 
(Ay — 6)? (4,0)? ‘ p 2s 0) 


* OSTROWSKI, A. M.: A method of speeding up iterations with super-linear con- 
vergence. J. Mathematics and Mechanics 7, 117—120 (1958) 
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Dividing this by (62) we obtain finally 

m } Pu 
A,.— 0 = (Gi; 

6 £2 #H=1 
( 3) (A,—0)3 = Dy 


v=1 Oy 


=O 


provided that the denominator on the right-hand side does not vanish. 


Observe that the limit on the right-hand side, valid for A)—>o, is in general 
different from the limit y in (46), valid for k-> oo. 


21. If, in the above iteration, the sequence A, converges to an eigenvalue o, 
we can obtain also a corresponding eigenvector as the limit of the vector sequence 
&,/|&|. Indeed, it follows from (30), (42) and (47) that if y,,..., y,, are the 
components of 7 corresponding to the eigenvalue o, then the sequence &,/| é,| 
tends to a vector €. The components of ¢ which correspond to the eigenvalues 
distinct from o vanish, while the components with indices 7,, ..., 7; are given by 


Wie Vig 
Lee VR e EH 


After a certain number of eigenvalues and corresponding independent eigen- 
vectors have been obtained by the above process, the method can be applied 
in such a way that in any case it gives either a new eigenvalue or at least a new 
eigenvector linearly independent of the eigenvectors already found. Suppose, 
indeed, that a complete linearly independent set of unimodular eigenvectors 


already found is given by ¢,,...,¢,;. Then choose a starting vector 7 orthogonal 
to all ¢,, 7.e., such that 
(64) (Gm) 0 i e)is 


Then I say that all vectors €&, deduced by the above process from 7 are ortho- 
gonal to each of the €,. Indeed, introducing normal coordinates as in § 2, we make 
f,,.--,&, the coordinate unit vectors corresponding to the indices 1,...,/. But 
then the orthogonality relation (64) reduces to 


Vie a 
and by (30) it follows that for each k 


that is, €, is orthogonal to each £,. Therefore the limiting value of &,/|&,| is also 
orthogonal to all ¢, and is independent of them. Thus our procedure started in 
this way gives us a new eigenvector independent of those already obtained. 

It is hardly necessary to mention that after » —1 linearly independent eigen- 
vectors have been obtained, the m-th eigenvector is determined by the n—1 
relation (64), and the corresponding eigenvalue is obtained at once as the value 
of the Rayleigh quotient. 

Further, it is easily seen to be sufficient that the relations (64) are satisfied 
with a certain approximation in order to be sure that the new eigenvector will 
not be a linear combination of the eigenvectors &,,..., &,. 
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22. In applying the iteration rule (28), (29), the largest amount of computa- 
tion is involved in the solution of the linear system with the matrix A —A,E, 
since this solution implies each time in the general case a multiple of n? elementary 
operations. It can be worthwhile therefore to reduce A by a preliminary trans- 
formation to a tridiagonal form, since for a system with a tridiagonal matrix 
the number of elementary operations implied in the solution of the corresponding 
linear system is only a multiple of n?, and the tridiagonal form is preserved in 
all the matrices A —A,E. On the other hand, as has been pointed out by 
W. GIVENS, the reduction of an Hermitian matrix to a tridiagonal form can be 
carried out by a finite number of “binary rotations’ *. 


However, the reduction to tridiagonal form would have to be carried out 
from the beginning with the greatest precision that may be required throughout 
the computation. Further, in this way we obtain directly only the approximations 
to the eigenvalues of the original matrix, since it probably would not be worthwhile 
to obtain from the eigenvectors of the transformed matrix those of the original 
matrix by applying again the whole sequence of binary rotations. And if we then 
obtain the corresponding eigenvectors of the original matrix A by applying the 
rule (5), we still have to solve the corresponding system of linear equations with 
an ill-conditioned general matrix. 


23. The eigenproblems in the case of a finite number of degrees of freedom 
usually present themselves not in the form 


(65) AE=HE 
but in the more general form 
(66) A&=ABE. 


In the form (66) B is a positive definite symmetric matrix while A is, as usually 
assumed, a real symmetric matrix. 


It is easy to deduce from our results the corresponding rules and theorems 
in the more general case (66). Observe that a positive symmetric matrix B can 
be written as a square of a symmetric positive matrix 4 


(67) aie. 


Introducing this expression of B in (66) and multiplying on the left by 47, 
we obtain 


AAA (Ae) 2 A (are 
therefore, putting 


(68) ATA AM A eS aes, 
(69) A,§, = Aé,, 


we see that any eigenvalue 4 of (66) is also an eigenvalue of (69) and the cor- 
responding eigenvector is obtained by (68) **. 


* A finite reduction to a tridiagonal form was contained implicitly in some previous 
papers by C. Lanczos but was not worked out explicitly. 


** Cf. CHERUBINO: Su un’equazione della teoria delle vibrazioni. Boll. Un. mat 
Ital., Ser. III, 11, 133—136 (1956) 
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If for an arbitrary vector £, the Rayleigh quotient is now formed relative 


to the matrix 4,, it can be written in terms of §=41&,, A and B by use of (69): 


Sitter Ae Cane) 
EE Sap e (eye 


On the other hand if we derive from &, a vector &()) by 
(A, — AE) fy =&,, 


this is reduced by (69) to 
(A — A, B) EM = BE. 


We see that our iteration rule (28), (29) in the case of (66) is transformed into 
the rule 


(70) &*) — (A — 4, B)7 Be, 
== Qa (Se) 
a tii 05 (Ey) 


and our main result (46) is therefore valid in the case (66) too*%. 
24. Examples**. 
ay le mattix. A == ( : A has 2 and 7 as eigenvalues. We begin with the vector 


eo (4, 0}, Obtain by (28) A, —=2,-thence compute’ s;)-Aj, .:.,'6,; 75, and finally 
apply to A,, A,, A; the improvement contained in the formula (61). The results 
are given in the following table: 


4 ii | ff) | a 
| | 
O 3 1 O 
1 2.076923 15 =| 
D 2.0419073413568826 | 2.079803921568627451 eal 
303. 2.05277 | 1.99836 | 1.9762747097570 a 
The convergence is particularly fast, since we have in this case y = 1G ae = 
1 ene 
jit Panay Iti O (pg 4.2) tem 0769; 
Sr The values of the quotients (i, 0)” ( O74 5-2); aoe are 0.0769; 


5 
0.0419; 0.0399; 0.05054. 


2. The matrix 


eas 
(2 5 6 | 
3 6 10 


has o,=0,=1 as a double eigenvalue, and o; =15 as the third eigenvalue. Start- 
ing with the vector &=(1, 0, 0), we obtain the values of Ag, &, Ay, Sa, Aa, 53, As 


* The rule (70), (71) is discussed directly in CRANDALL’s paper cited above. 
*« For the computation of these examples I am greatly indebted to Mrs. Iba 
Ruopes, National Bureau of Standards, Washington D.C. 
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and derive again from A, A,, A; the value of A; by (61). The results are contained 
in the table: 


i | ® EVES 
ae oe Fae: 
1 | 1.0°6369426751592356687898089171974 | —12 pam 
2 | 1.081320194335880903168602516 — 13.0°4901 96078431372 2 | 3 
B33" 1.0811739 | 0.93183 3712.9 53776 2h 33 
The value of 28—% and those of “##2—% (k =O, 1, 2) are in this case 0.001 448, 
A,— 6 (A,— 9)? 
0.006369, 0.005109, 0.005 1017 while y = 58 is 0.005 10204. 


American University, Washington D.C. 
and 
University of Basle, Switzerland 


(Received September 29, 1958) 


A Short Proof of a Boundedness Theorem for Linear 
Differential Systems with Periodic Coefficients 


JEK HALE 


Communicated by L. CESARI 


Introduction. In the last few years a number of boundedness theorems for homo- 
geneous, linear, differential systems with. periodic coefficients and containing a 
small parameter have been proved by L. CEsart, the author and others, by using a 
method of successive approximations first devised by CEsArI. More recently, 
M. GoLomB [6] has proved similar theorems by using another method of succes- 
sive approximations devised by WEBER (see [2], in particular § 4, no. 5, for 
general references). The aim of the present paper is to show that an elementary 
argument leads to a new general theorem of boundedness which comprehends 
as a particular case some (though not all) of the statements previously proved 
by the method of CesariI and by the method of WEBER. 


Consider the system *: 


2 
u; + A;u, = 2 [Pint A) 4, + y(t, A) up] + Vig (t,A) us, 7 =1,2, 
(1) 


: 
uz = >) [Dyy(t, A) my, + Poy (t, A) ta] + Bea (t, A) us, 


> 
ll 
ran 


where / is a real parameter, 3 1s a scalar, u, is an m-dimensional vector, O<m<n, 
Uy 1s an n —m-dimensional vector, n=m, and each ®;,(t, 4), %,,(é, A) is a matrix 
whose dimension is determined by (1) and the elements of these matrices are 
real-valued functions of the real variable ¢, periodic in ¢ of period T =2z/a, 
L-integrable in [0, 7], and are continuous functions of A at A=0 for almost 
all ¢ in [0, 7]. If p(é, A) is any element of these matrices, then we assume that 
p(t, 0) =O and there exists a function q(t), L-integrable in [0, T], and a A, >0 
such that | (¢, 4)|<g/(t) almost everywhere in [0, T], OS|A|SA). The class of 
functions p(t, A) which satisfy these properties will be called A(7). Furthermore, 
A, = diag (oj, -.., 0%,), Ay = diag (63,4, .--; 04) where o,, ..., 0, are positive real 
numbers with o;+0,=-7@, 20;=70, 0,7, ES Hao 6 ae Oe = 
OMG ES aac. 

We shall say that a function f(t) is even with respect to y if f(t+y) =f(—2) 
and that it is odd with respect to y if f(t +yv) =—f(—?). 


* Throughout the present paper, two functions are said to be equal if they are 
equal except for a set of Lebesgue measure zero. 
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As it is immediately seen, these two concepts coincide with the usual evenness 
and oddness when the translatation ¢=t,—4y is made. For real functions /(?) 
periodic of period T=2x/m, L-integrable in (0, 7], with Fourier coefficients 
Ay, 4, 0,,k =1, 2,..., the following cases are of interest: f(t) is even with respect 
to (a) y=o if f®=f(—4), (b) y=sT, if 4,-1=%2,=9, RNC eee 
(c) y=7/4k if {() =a(cosk@t-+ sinkat); f(Z) is odd with respect to (a’) y =0, 
if f(t) =—f(—d, (b') y=T/2 if ag, =02,1=0, R=1, 2,.-., (ce) y= Taka 
}(t) =acos2kwt+bsin4kwt. We prove the following: 


Theorem 1. If the matrices D;,,'%j, of system (1) are such that ®;,(—1t) = 
(—1)*#7O,,¢+y), Ya(—f) =(—1)7* E+) for some real number y, then 
there exists a A, > 0 such that all of the absolutely continuous (AC) solutions of (A) 
are bounded for — 0 <t<+ oo, OS|A| =i. 

In case all of the coefficients in the last equation of system (1) are zero and 
the coefficients of v3 are zero in the first m equations of (1), then Theorem 1 is 
a boundedness theorem for systems of second order differential equations. For 
this case, Theorem1 for y=0, m=mn and all the matrices Y.=0 was first 
proved by L. Cesar [3] and was later extended by R. A. GAMBILL [5] to the 
case where %,-+-0. For y=0, m=, Theorem 1 was first proved by the author 
[8] and later, for all %,—0, by M. Gotoms [6]. Then the elements of the 
matrices ®,,, Y%, are either even functions of ¢ or odd functions of ¢. 

The elementary argument used in this paper for the proof of Theorem 1 
seems to be new, though some elements of it can be traced to A. Lyapunov [10, 
p. 409] for systems of order m and to L. CEsAr1 & J. K. HALE [4] for systems 
of order m=2. On the other hand, by using the method of Cesar, the author 
[9] has been able to obtain a boundedness theorem which includes Theorem 1 for 
systems of type (1) where more than one characteristic root for A=0O may be 
zero. The emphasis in this paper is to show that many results concerning the 
qualitative behavior of the solutions of systems of type (1) can be obtained 
without using successive approximations. 


Proof. In this section we prove Theorem1, but we begin by discussing the 
properties of a system of differential equations which is slightly more general 
than system (1). Consider 

2 = AZ + Cit, 2) 2 (2) 


where / is a real parameter, z=(%, ..., 22,43), A is a (2n +4) X(2m4+4) matrix 
whose elements are real constants and C(t, A) is a (2n+4)x(2m-+14) matrix 
whose elements belong to the class A(T). Furthermore, we shall suppose that 


the characteristic roots v,, k=1, 2,..., 2n +4, of A are Vg j;-1=16;, Vzj= BOR 
Yenti=9, where o,, ..., 0, are positive, real numbers such that UV; == UV, (mod w1), 
G2zh, 4G, RA, 2;...,20-+4. 
Let 
d = min |im w — (v;— v,)|, 


Pome Qi noes 2 et. 
m=0,1,..., |7—k| +-m>0 


and let ¢ be any real number such that 0< é€<¥, and consider, in the complex 
plane, circles C,,...,C,,4, with radii ed such that C,; has center at oj7 =1, 
2,-.-,2M+41. If wy, ..., Mens are any complex numbers such that pp ECy, 
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in 2 hie. ree Aevithen 


fom ey ee |imw — (u; — u,)| 24d — 2ed =d(1 —2e)>0. (3) 
m=0,1, 2,...,|7—k| +m>0 

Let t;=1;(A), 7=1, 2,...,2%+1, denote the characteristic exponents of 
system (2). The numbers 1,, of course, are only determined up to a multiple 
of m7, but we may assume without loss of generality that t,;_,(0) =7o,, T2; (0) 
=—10;,]7=1, 2,...,%, Te,41(0)=0. Furthermore, exactly as in [4], it can be 
shown that the above conditions on the matrix C(t, 4) imply that these charac- 
teristic exponents for system (2) are continuous functions of 2 at 4=0 and 
there exists a Aj>O such that for OS|A|SA, t.;_1(A) =72;(A), 7 =1, 2, ..., 0, 
and 7,1 is real. Also, since t,(A) is continuous at A=0, there exists a Aj >0 

such that 1;(a)€C,; for OS|A|</o, which implies from (3) that 


t; (A) + 7, (A) (moda), Poe Ke Agha 2D, pangs es | (4) 


Therefore, we have the following: 


Theorem 2. For 0S|A|SA,, Ap=min (Ap, Ag, Ao), the characteristic exponents 
T;(A), 7=1, 2,...,2m-+41, of system (2) have the following properties: 


(i) T9271 (A) SF Tz; (A) , ] ae 1, 2, +++), Tant1 (A) iF Tan+1(A) , 

(li) 7; (A) 32 t, (A) (mod wz), ROK, IR =e NG lee ated 

If Bch, $= 1. Qi. Oh and OS|A| <A/,, then, it follows from Theo- 
rem 2 that the f; are distinct complex numbers. It is known that, when all of 
the characteristic multipliers 6; are distinct, there is a fundamental system of 
AC solutions of (2) of the form 


eget > k= 1, Qin, 20th fe, 2.0.2, 2H rp, (5) 


where each #,,(t) =—;,(¢+ 7). 

Theorem 3. Let z(t) be any solution of (2); let P be a real constant nonsingular 
matrix which is independent of A; and let (*) denote the differential equation for 
w(t) = Pz(—t) which is obtained from (2) by replacing t by —t and applying the 
transformation P. If there exists a real constant nonsingular matrix Q which 1s 
independent of A and such that w(t)=Q2z(t) also satisfies (*), then all of the AC 
solutions of (2) are bounded for — o<t<+oo, 0OS|A|<d,, where 2, is defined 
above. 

Proof. For 0S|4|S/,, Z() ={%, (0,7, R=1, 2, ...,2%+1} defined by (5) 
represents a fundamental system of AC solutions of (2). It is clear that PZ (—?) 
is a fundamental system of AC solutions of (*). If z(¢) is any solution of (2), 
then there exist a constant vector a such that 


z(t) =Za. 


We choose z(é) so that each of the components of the vector a are different from 
zero. On the other hand, since Qz(t) satisfies (*), there exists a constant vector 


6 such that 
Q2(t) = PZ(—2) 0. 
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But, the previous formula together with this one implies that 
OZ YEE Za 1) e- 


However, from Lemma (2.ii) of [7], this formula implies that there exists, 
for every 7=1,2,...,2%+1, an /(j) such that —1t,=1, (mod w1); or, if R(a) 
and 3(a) denote the real and imaginary parts of a, respectively, 

— R(z,;) = R(t,) (moda), 

— 93 (t;) =S(t) (moda). 
From Theorem 2, t, must be the complex conjugate of t; and, thus, R(t;) =0, 
j=1,2,...,2n+41. Finally, from (5), all of the AC solutions of (2) are bounded 
for — 0<t< oo, 0S|A|<A, and the theorem is proved. 

If the vectors 4, , %, in (4) are given by #,= (u41,.--, Uym), Uo= (Wor, «++, Ua, n—m)> 
then the transformation of variables 

2 == Uy. Zn4j = M1; 1 SAD ies ed 
hs NO et 1.=152)5.. pH Sm; (6) 
Zon +1 as 


leads to the equivalent system of 2” +1 first order equations 


2 =Az+Ctt, d)z, (7) 
where 
A v4 Va By, ) G ze es oe 
By, OF aan @ wy (8) 
— A, OF n—m 
Byg = (L,,0n1), Iepre= Oa —A, 4 
O07, m 0;, m—™M 


O;% is the 7 xk zero matrix, I, is the »Xm unit matrix, O=(@,, 
where the matrices ®,, %, are defined in (1). If we let w(t) = Pz 
the matrix P is defined by the transformation 


eT, 
S 

= 

a 

iv) 

ie) 

a) 


w,(t) =2,(—12), O45 (t) = —2y4(— 8), pied |. Noo 
w(t) = = oe t) , y+; (¢) = 2a t) , ] Sd Pia aks ea, 
Wey4i(t) = — Zonti(— 4), 


then the differential equation (*) in Theorem 3 is given by 


w =Aw-+C*i(t, a), (9) 
where 
*k OS OF 
C -( fee ) (10) 


S 
* 
Ss 
= 
I 


[(—1)**7@,,(—#, a), 7 =1,2,3, R= so, 
P(t, A) = [(— 1) (6A), 7, k =1, 2, oe 
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From the remark preceding Theorem 1, we may assume y=0. If the assumptions 
of Theorem 1 are satisfied for y=0, then, ®*(¢,2)=@(¢, 4), Y*(t, 2) = W(t, a). 
Therefore, Theorem 1 follows from Theorem 3 with Q =/,,,4. 
The most natural question to ask is whether it is possible to find matrices 
P, Q which are essentially different from the ones constructed for the proof of 
Theorem 1 and, therefore, obtain a boundedness theorem more general than 
Theorem 1. To investigate this question, it is clear that Q may be taken to be 
Z,,,4, and it will be sufficient to find a non-singular matrix P which is independent 
of A such that 
P[A + C*(t, A)] P1=A + Ct, a) (11) 


for all ¢ and A, where C* is defined by (9) and (10) and C is defined by (8). It is 
convenient to let P=(P,), 7, =1, 2,3 where each P, is a matrix and B,, By 
are m xn and Fy, is a scalar. Since (11) must be satisfied for all 2, this implies 
Cate A ornate, = Lge Ae = ar{ Ay, A,))b4, £.diagonal, 
P,, arbitrary and all other P, are zero matrices. Furthermore, if each of the 
equations in (1) contains at least one nonzero periodic coefficient, then it is easy 
toxshow that Fy. is the-zero matrix. Jf we let Bj —diag(p,, .....p,). Rs=Pasa 
NG = Liars Wir Ay BA; Dyce I) pei Wig Iu HAs 2ye ney Met 1) a then 


J 
the relation PC*(t¢, 4) =C(t, A) P leads to the system of equations 


P; Pjn(— 4,4) =4 Ob. G2 (6 A), Gees erent oe Ane = Na Dene eI 
DP; Vjx(—t, A) —— =n Dr Vel, A), 4, R lee An iy (LP +1, 


for the undetermined quantities f,,...,,,,.. lf we let —*=% im the first 
equation of (12), then we have 


, AYP 
yl A) = tO gpa %, 2) = (BY gala) 


(12) 


Pj 
and, therefore, (p,/p;)?=1, 7, =1, 2,...,% +4. By a change of the scale of the 
independent variable in (1), one can assume that each #;,7=1, 2,...,2u +1, 


is either +1 or —1. But, with each ; either +1 or —1, it is clear that there 
exists an elementary transformation which will make conditions (12) equivalent 
to the conditions in Theorem 1. Therefore, in this sense, Theorem 1 is the most 
general theorem for this class of functions which can be obtained by using 
Theorem 3. 


The majority of this work was completed while the author was at Remington 
Rand Univac, St. Paul, Minnesota. 
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Finite Integrity Bases for Five or Fewer 
Symmetric 3x3, Matrices 


A. J. M. SPENCER & R.S. RIVLIN 


1. Introduction 


In a previous paper [/], finite integrity bases for five or fewer symmetric 
3X3 matrices, under the orthogonal transformation group, have been derived. 
In the present paper, it will be shown that the integrity bases there derived are 
redundant in the sense that certain of their elements can be expressed as poly- 
nomials in the remaining elements, and integrity bases containing fewer elements 
will thus be obtained. The results of the previous paper [7], which form the 
starting point of the present paper, are essentially consequences of the Hamilton- 
Cayley theorem and the procedures adopted in the present paper are similarly based. 

Finite integrity bases for five or fewer symmetric 3 x3 matrices under the 
general linear transformation group have been given by Topp [2, 3] and Turn- 
BULL [4]. Since, by PEANo’s theorem, a finite integrity basis for any number 
of symmetric matrices can be derived by polarization from an integrity basis 
for five symmetric 3 x3 matrices together with the alternating invariant of six 
symmetric 3 x3 matrices, Topp’s result can be used, as he pointed out, to con- 
struct a finite integrity basis for any number of symmetric 3 x3 matrices under 
the full linear group. In particular it can be used to construct a finite integrity 
basis for six symmetric 3 x3 matrices under the full linear group. By replacing 
one of these matrices by the unit matrix, a finite integrity basis for five symmetric 
3 X3 matrices under the orthogonal group can be obtained in accordance with 
the Adjunction theorem (see, for example, WEYL [5]). The elements of the finite 
integrity basis thus obtained from Topp’s results may, with considerable algebraic 
manipulation, be expressed as polynomials in traces of matrix products and a 
number of these traces which form an integrity basis may be selected by inspec- 
tion. Certain of these may be expressed as polynomials in the remainder by 
employing the Hamilton-Cayley theorem, or results obtained therefrom by 
appropriate substitutions, and the end results of the present paper can be 
obtained. However, since the labour involved in obtaining these results from 
those of Topp is comparable with that of deriving them directly, we have adopted 
the latter course in the present paper. 


2. Finite integrity bases for five or fewer matrices 


In this section we give the integrity bases for one, two, three, four and five 
symmetric 3 X3 matrices, under the orthogonal transformation group*, which 


* It is immaterial whether we consider the full or the proper orthogonal group. 
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were obtained in the previous paper [/]. Some of the results, e.g. those for 
a single matrix and for two matrices are, of course, well-known. 


Single matrix a. ca) tere, tra’. (2.1) 
Two matrices a, b. tr a, tr a”, tra’, 
tr b, tt b?: melas (2.2) 


trab, trab?, tra?b, tra*b?. 


Three matrices a,b,c. The integrity bases for each of the three pairs of 
matrices which can be selected from a, b, c, together with the invariants 


trabe, trabe? and trab?c? (25) 


and the invariants which can be formed from (2.3) by permutations of a, b, €. 


Four matrices a, b,.c,d. The integrity bases for each of the four sets of 
three matrices which can be selected from a, b, c, d, together with the invariants 


trabed, trabcd?, trabe?d?, trbcada@ and trbacd@& (2.4) 


and the invariants which can be obtained from these by permutations of a, DB, e, d. 


Five matrices a, b,c, d,e. The integrity bases for each of the five sets of 
four matrices which can be selected from a, b, ¢, d, e, together with the invariants 


trabede, trabcde?, traebcde?, trabecde?’, 


De 
trabcede2 and trabcde? eee 


and the invariants obtained from these by permutations of a, b, e, d, e. 


3. Some theorems concerning traces of matrix products 


In § 4 and the following sections, we shall show that certain of the elements 
of the integrity bases for three, four and five matrices given in § 2 can be expressed 
as polynomials in the remaining elements. In this way we shall derive integrity 
bases for three, four and five matrices containing fewer elements than those given 
in § 2. In order to do this we shall require certain results which are either stated 
or derived in a previous paper [/]. These are given in the present section. 


Firstly, we require the following two lemmas which result immediately from 
the definition of matrix multiplication. 


Lemma 1. The trace of a matrix product formed from 3 <3 matrices is unaltered 
by cyclic permutation of the factors in the product. 


. Lemma 2, The trace of a matrix product formed from symmetric 3X3 matrices 
vs unaltered if the order of the factors in the product 1s reversed. 


If J, and J, are two scalar invariants of any number of 3 X3 matrices of 
equal partial degrees in each of the matrices and I,—J, is expressible as a poly- 
nomial in invariants of lower or equal partial degrees in each of the matrices 
and of lower total degree, we say J, and J, are equivalent and write 


1, = Ty. 
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We also say I,—TI, is reducible and write 
I,—I,=0. 


We require the following relations which, in the previous paper [7], are seen 
to result from the Hamilton-Cayley theorem. Let a, b, ¢ be three 3 x3 matrices 
and let y and 2 be any 3 x3 matrices (in particular cases I). Then, we have 


tryabez+trybcaz+ trycabz+ 

“4 
+trybacz+tryacbz+4trycbaz=0, 6-1) 

unless y=2=T, and 
tryaba’z+try@baz=0, (3.2) 

even ity —< =f, 
Also, if # (+=) is a 3 X3 matrix, we have 

try aa b?z=0, (3.3) 


unless ai 2.— 1. 


4. Reduction of the integrity basis for three matrices 
We take as our starting point the finite integrity basis for the three sym- 
metric 3x3 matrices a, b,c, given in §2, and consider equivalences of the 
invariants (2.3) and the invariants formed from them by permutations of a, b, e. 
From lemmas 1 and 2, it is immediately obvious that trabe is unaltered 
by permutation of a, b, e. 
In a similar manner it follows from Lemmas 1 and 2 that 


tra b?e? = trac?b? 
and 
trabc?=trbaec’. 


We thus see that an integrity basis for the three symmetric 3 x3 matrices 
a, b,c under the full or proper orthogonal group is formed by the integrity 
bases, given by (2.2), for the three sets of two matrices which can be selected 
from a, b, ec, together with the invariants 


trabic.. trabdc.. ptrbca* s treeb:,. trab*e, un 
4.1 
trbc?a? and trea’®b?. 


5. Reduction of the integrity basis for four matrices 

The elements of a finite integrity basis for four symmetric 3 3 matrices 
a, b, c,d, which involve all four of the matrices, are given by (2.4) and the 
invariants formed from these by permutations of a, b, e,d. We shall consider 
in turn each invariant listed in (2.4) together with the invariants formed from 
it by permutations of a, b, c, d, with a view to determining which of these may 
be omitted from the integrity basis in the light of the results given in § 3. 

trabed. From Lemmas 1 and 2 we see that each of the invariants formed by 
permutation of a, b,¢,d in trabed is equal to either trabed, trabde or 
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tracbd. From (3.1), replacing a, b, ¢ by b, ¢, d respectively and taking y =a, 
2 =I, we obtain 
2(trabed + trabde + tracbd)=0. (5.4) 


Thus, we can omit one of the three invariants on the left-hand side (say tr acbd) 
from the integrity basis and, of the terms tr abed and the terms formed from 
them by permutation of a, b, ¢, d, include only 


trabed and trabde. (5.2) 


tr abed?. From lemmas 1 and 2, we see that each of the invariants formed 
by permutation of a,b,c in trabed? is equal to either tr abed?, tracbd? 
or tr bacd?. From (3.1), taking y=I, e=d?, we have 


2(trabcd*+ tracbd’+trbacd’)=0. (5.3) 


Thus, we may omit one of the invariants on the left-hand side of (5.3), say 
tr bacd?, from the integrity basis and, of the terms tr abed? and invariants 
formed from it by permutations of a, b, €, include only 


trabcd* and tracbd’. (5.4) 


Applying similar arguments to the invariants formed by permutations of d, a, b 
in tr dabe?, of c,d, a in tr cdab? and of b, c, d in tr beda?, we see that of the 
invariants obtained from tr abed? by permutations of a, b, c,d, we may omit 
from the integrity basis all except 


trabed?, tracbd?, trdabec?, trdbacec?, 


tredab?, treadb?, trbeda, trbdca’. 6.5) 


tr abe?d?. Now, we consider the invariant tr abe?d? and invariants formed 
from it by permutations of a, b, e, d. 


Replacing ¢ and d by ec? and d? in (5.1), we obtain 
trabe’d+trab@e?+traecvow=o. (5.6) 
Now, taking y=a, =I in (3.3) and replacing a, a, b by c, b, d respectively 
we have 
traebd=0. (5.7) 
From (5.6) and (5.7), we obtain 


traberd@= —trabede?. (5.8) 
Employing Lemmas 1 and 2, we obtain 
treb@e’=trbac’d’. (5.9) 


From relations of the types (5.8) and (5.9), we see that of the invariants 
formed by permutations of a, b, ¢, d in tr abc?’ d?, we may omit from the integrity 
basis all except 

trabe’'d, tracb?d, trad b2e?, 


trbca@’d@, trbda@ec?,  treda’b?. bond) 
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tr bacda@ and tr beada?®. Finally, we consider the invariants tr bacda? 


and tr beada? and the invariants obtained from them by permutations of 
a, b, e, d. 


Taking y =I and 2 =da? in (3.1), we obtain 


tr(abc+acb+bac+bea+cab+cbad@=o. 
Since trabcda?=0 and tracbda?=0, we have 
tr(bacda@ + beada+cabda?+cbada’)=0. (5.11) 
Now, from Lemmas 1 and 2, we have 
trbcada@=trad@bec=trebarda. (ee) 
From (3.2), replacing b by d and taking y=cb, =I, we obtain 
treba@da+trebad@’=o0. (5.13) 
From (5.12) and (5.13), we obtain 
trbcada@’= — trebada’. (5.14) 
Employing the relation (5.14) in (5.11), we obtain 
trbacda= —treabda’. (5205) 
Again employing Lemma 1 and the relation (3.2), we obtain 
trbacd@=trced@ba=—tredaba. (5.16) 


From (5.16), we see that we may omit from the integrity basis tr beada?® and 
invariants formed from it by permuting a, b, ¢, d. 


From (5.14) and (5.16), we have 
trbacda@= —trbadea@. (5.17) 


From (5.15) and (5.17) we see that we may omit from the integrity basis all but 
one of the set of invariants consisting of tr bacda® and the invariants formed 
by permuting Jb, c, d in this. 

Thus, of the invariants tr bacda? and tr beada’ and the invariants obtained 
from these by permutations of a, b, c,d, we may omit from the integrity basis 
all except 


trbacda, trebdab*, trdcabc® and tradbecd. (5.18) 


Collecting the results obtained above, we see that an integrity basis for the 
four 3X3 symmetric matrices a,b, c,d under the full or proper orthogonal 
groups is formed by the integrity bases for the four sets of three matrices which 
can be selected from a, b, ec, d, together with 


(i) trabed, trabdc; 


(ii) trabed, tracbd?, trdab c?, trd bac?, 
tredab?, tread b?, trbeda’, trbdca’; 
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(iii) trabe’d, tracb?d?’, trad b*c?, 
trbca*d’, trbda’c?, tred a? b?; 
(iv) trbacda, trebdab’, trdcabe’ and trad bed’. 


6. Reduction of the integrity basis for five matrices 
In this section we shall discuss the equivalence of the elements, given in § 2, 
of the finite integrity basis for the five symmetric 3 x3 matrices a, b, €, d, e, 
which involve all five of the matrices. The procedure will be similar to that 
adopted in § 4 in discussing the case of four matrices. 
trabede. We consider first the invariants formed from trabede by 
permutations of a, b,c,d,e. Taking y=de and z=I in (3.1), and employing 
the notation 
NMabc=abc+acb+bea+bac+cab-+cba, (6.1) 
we obtain 
tr(de>\abc) =0. (6.2) 


Noting that >'abe is a symmetric matrix, we have, with Lemmas 1 and 2, 


tr(de abc) =tr(ed > abc) =tr{e(. abe) d} 


=tr{d().abe)e}=0. (6:3) 
We can form ten relations of the type (6.2). These are 
tr(de>\abc)=0, tr(cee>abd)=0, tr(ed> abe)=0, 
tr(be>,acd)=0, tr(bd>\ace)=0, tr(be>dade)=0, 
tr(ae>| bed) =0, tr(ad>}bce)=0, tr(ac> bde) =0, es 
tr(ab > cde) =0 


Not all of the relations (6.3) are independent, for it may be readily verified by 
introducing relations of the type (6.1) that 


tr(ac > bde) + tr(ad > bce) + tr(ae> bed) 

= tr(be Dade) + tr(bd S.ace) + tr(be acd), 
tr(ab >’ cde) + tr(ad >) bce) + tr(ae > bed) 

=tr(cb Diade) + tr(ed S abe) + tr(ce abd), 
tr(ab > cde) + tr(ac > bde) + tr(ae > bed) 

= tr(db Dace) + tr(de Yabe) +tr(ded abo), 
tr(ab > cde) + tr(ac > bde) + triad Y bee) 

= tr(eb Diacd) + tr(ec Yabd) +tr(ed Yabe). 


(6.5) 


From equations (6.5), we can derive a number of similar equations. For 
example, subtracting the second equation from the first, we obtain 


tr(ba > ede) + tr(bd ace) + tr(be Yaed) 


= tr(ca > bde) + tr(ed abe) +tr(ce Yaba). ioe) 
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The four equations (6.5), regarded as equations in which the unknowns are 
tr(ab>\ ede), tr(ac> bde), tr(ad> bce) and tr(ae> bed), are independent and 
may therefore be used to determine these quantities in terms of tr(bc>\ade), 
tr(bd> ace), tr(be> acd), tr(cd>\abe), tr(ce>\abd) and tr(dey abe). Thus, 
the last four of the relations (6.4) may be regarded as consequences of the first six. 

These first six relations may be re-written, using Lemmas 1 and 2, as 
trabcde-+trabced+trabdec+trabedc+tracbde+tracbed=0, 
trabced+trabdce-+trabdec-+trabecd+tracebd-+tradbce=o, 
trabcde+trabdce+trabecd+trabedc-+tracdbe+tradcbe=0, 
trabecd+trabedc+tracbed+tracdbe-+tracebd-+tradcbe=0, 
trabdce+trabdec-+tracbde+tracdbe+tracebd+tradbce=0, 
trabcde+trabced-+tracbde+tracbed-+tradbce+tradcbe=o. 


(6.7) 


Regarding these equations as a set of six simultaneous equations in which the 
unknowns are trabced, trabdce, trabedec, tracbde, tracebd and 
tr adcbe, we can show, by considering the rank of the matrix of the coefficients 
of the unknowns, that they can be used to determine the six invariants in terms 
of trabede, trabdec, trabecd, tracdbe, tracbed and tradbce and 
invariants involving four or fewer of the five matrices. Thus, of the invariants 
formed by permutations of a, b,c, d,e in trabcde, we may omit from the 
integrity basis all except 
trabcde, trabdec, trabecd, 


(6.8) 
tracdbe, tracbed, tradbece. 


The set of invariants (6.8) are all formed from tr abede by even permutations 
of the matrices a, b, c,d, e and it is easily shown that any invariant formed 
from tr abede by even permutation of a, b, ec, d, e may, by means of Lemmas 1 
and 2, be expressed in one or other of the forms (6.8). The choice of invariants 
to be retained in the integrity basis is by no means unique; for example, in place 
of the set (6.8), we might have chosen to retain the invariants formed by odd 
permutations of a, b, c,d, e in trabede, and other choices are also possible. 

trabcde?. We next consider the invariants formed from trabede? by 
permutations of a, b,c,d,e. 

Replacing e by e? in the above discussion, we see that each of the invariants 
formed from trabecde? by permutations of a, b, c,d can be expressed as a 
pa trabcde*, trabde’c, trabe’cd, 


6.9 
tracdbe?, tracbe’d, tradbce?’. io) 


Thus, noting that trabde?c=treabde’, trabe’cd=tredabe’ and 

tracbe?d —tr dacbe?, we see that of the invariants formed from tr abcde? 

by permutations of a, b, ce, d, e, we may omit from the integrity basis all except 
trabede?, treabde?, tredabe?, . 


6.10 
tracdbe?, trdacbe’, tradbce’, ( ) 


and invariants formed from these by permuting a, b, ¢, d, e cyclically. 
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trabcede?, tr abecde?, traebcde*. Now, consider the invariants formed 
by permutations of a, b, c,d, e in trabcede?, trabecde® and traebcde’. 

Replacing a and b in (3.2) by e and aed respectively and taking y—b, 
z=I, we obtain, with Lemma 1, 


trbeacde’=— trbeacde = — tracde be’. (6.11) 
It follows from (6.11) and relations obtained from it by permuting a, b,¢,d 
that, of the invariants formed by permutations of a, b, c,d in trabcede?, 


trabecde? and tr aebcde?, we may omit from the integrity basis all except 
those formed by permutations of a, b,c,d in trabecde?® and trabcede’. 


Replacing a and b in (3.2) by e and ed respectively and taking y=ba, 
z =I, we obtain, with Lemma 1, 


trbaecde?= — trbae’cde = — tredebae?’. (6.12) 

Using Lemmas 1 and 2, we can obtain 
iredebae? =tre’cdeba =trabedce’. (6.13) 
From (6.12) and (6.13) and relations obtained from these by permuting a, b, ¢, d, 
we see that, of the twenty-four invariants obtained by permutations of a, b, ¢, d 


in tr abecde?, we may omit from the integrity basis all except six, which may 
be chosen as 


trabecde’?, trbaecde?, tracebde?’?, treaebde?, 


| 
trbceade*® and trebeade?. ot) 


Replacing a, b, cin (3.1) by a, e, d respectively and taking y = J and z = bce?, 
we obtain 
traed+ade+dae+dea+ead+eda)bece=0. (6.15) 


Since tr eadbce?=0 and tr edabce?=0, the relation (6.15) yields 
traedbce’+ tradebce?+ trdaebce?+ trdeabce=0. (6.16) 
By permutation of a, b, c,d in (6.11), we obtain 


traedbce’=— trdbceae’? and trdeabce?=~— trabcede. (6.17) 
From (6.16) and (6.17), we have 


tradebce’+ trdaebce®— trdbceae?— trabcede=0. (6.18) 
From (6.12) and (6.13), we have 


trabedce’?= — trbaecde?. (6.19) 
By permutations of a, b, e,d in (6.19), we obtain 
tradebce’=—trdaecbe® and trdaebce?= — tradecbe.. (6.20) 
Permuting 6 and © in (6.18) and employing the relations (6.20), we obtain 


—trdaebce’— tradebce— trdcbeae— tracbede=0. (6.21) 
The relations (6.18) and (6.21) yield 


trdbceae’+ trabcede®?+ trdebeae+ tracbede=0. (6.22) 
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From Lemmas 1 and 2 we have 
trabcede?—trdecbae®? and tracbede?=trdebcae®. (6.23) 


With (6.23), we may re-write (6.22) as 


tr{d(bece+ecb+cbe+ebcjae}=0. (6.24) 

Now replacing a, b,c in (3.1) by b,¢,e respectively and taking y=d and 
2 =ae?, we obtain 

tr{d(bcee+ecb+cbe+ebc+ceb+becjae}=0, (6.25) 


which with (6.24) yields 
trdcebae*? = — trdbecae’. (6.26) 


Employing relations obtained from (6.26) by permutations of a, b, c,d, we 
see from the result (6.14), that, of the invariants formed by permutations of 
a, b, c,d in trabecde?, we may omit from the integrity basis all except 


trabecde’, trbceade* and treaebde?. (6.27) 


Taking all possible permutations of a, b,e,d in the relations (6.18), we 
obtain, with the result (6.27), the relations 
trabcede’?+ trdbceae’= q,, 
tracbede’+ trdcbeae’= qy, 
trebaede’®+ trdbaece’®= gy, 
tradbece’+ tredbeae’=qy, 
trbacede’?+ trdacebe?=9,, 
trbadece’+ trcadebe?= qe, (6.28) 
tracdebe*?-+ trbcdeae’=_q,, 
tradcebe*+ trbdceae’= gq, 
treabede?-—+ trdabece’= gq, 
trabdece’?+ trebdeae= Op, 


trbcaede?+ trdcaebe? =, 


trbdaece®+ tredaebe= yp, 


where the q’s are linear combinations of the invariants (6.27). 
Not all of the relations (6.28) are independent. For example, adding the first 
and third of the relations (6.28), we have 


trabcede?+ trebaede?+ trdbceae’+ trdbaece’=9g,+ 93. (6.29) 
From Lemmas 1 and 2, we have 
trebaede=trede’cba=trabcede. (6.30) 
Thus, with a relation of the form (3.2), 
trabcede?+ trebaede*=trabc(ede®+ ede) =0. (6.31) 


50x 


44-4 A. J. M. SPENCER & R. S. RIVLIN: 


Introducing (6.31) into (6.29), we obtain 


trdbceae®+trdbaec@= q+ o3- (6.32) 

Again, from relations of the form (3.2), together with Lemmas 1 and 2, we have 
trdbceae’= — trdbcetae = — treaedbe=— trebdeae’, 

and (6.33) 
trdbaece= — trdbae’ce = — trecedba = — trabdece’. 


From (6.33) and (6.32), we have 
trebdeae?+ trabdece’= — Gy — 9s, (6.34) 


which expresses the same fact as the tenth of the relations (6.28). Thus, the 
tenth of the relations (6.28) can be derived from the first and third. In a somewhat 
similar manner we can derive the twelfth relation from the fourth and eigth; 
we can derive the ninth relation from the fifth and sixth and the eleventh relation 
from the second and seventh. Thus, the last four of the relations (6.28) are 
derivable from the first eight. 


By appropriate permutations of a, b, c,d in (6.33), we obtain 


trdcbeae’= — trbcdeae’, 
trabcede?=—trebaede?’, uae 
trdacebe?= — treadebe?, OS 
d 
se trbdceae?= —tredbeae?’. 


Substituting from (6.35) in the first eight of the relations (6.28), we obtain 


—trebaede’+itrdbceae’= 9, 
tracbede*’— trbcdeae’= yu, 
trebaede®+trdbaece’= q,;, 
tradbece®?+ tredbeae= Qy,, 
trbacede?— trcadebe?=q,, 38 
trbadece®?+ trceadebe= q,, 


tracdebe?+ trbcedeae=~gq. 


7a) 


tradcebe?— tredbeae?= Qz. 


Now, the twenty-four invariants obtained by permuting a, b,e,d in 
tr abcede?® are the twenty-four traces occurring on the left-hand sides of the 
relations (6.28). From (6.36) it is seen that each of the sixteen invariants occurring 
in the first eight of these relations is equivalent to some linear combination of 


trbcdeae’, trebaede®?, trcadebe? and tredbeae (6.37) 


and the three invariants (6.27). By means of relations obtained from (6.35) by 
appropriate permutations of a, b, c, d, we see that each of the traces appearing 
in the last four of the relations (6.28) is equivalent to the negative of one of 
those occurring in the first eight relations. We therefore have the result that 
each of the invariants obtained by permuting a,b,c,d in trabcede? is 
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equivalent to a linear combination of the invariants (6.37) and (6.27). Together 
with the result following from (6.11) and the result (6.27), we see that, of the 
invariants formed by permutations of a, b,c¢,d in trabcede?, trabecde? 
and tr aebcde?*, we may omit from the integrity basis all except the invariants 
(6.37) and (6.27). 

Consequently, of the invariants formed by permutations of a, b, ¢, d, e in 
trabcede®, trabecde? and traebede?, we may omit from the integrity 
basis all except the invariants (6.37), (6.27) and the invariants obtained by cyclic 
permutation of a, b, c,d, e in these. 

trabed’e*. Finally, we consider the invariants formed from tr abed?e? 
by permutations of a,b,c,d,e. 


Taking y=IJ and z=d?e? in (3.1), we obtain 
tr{(abe +acb + bca+bac+cba-+cab)de}=0. (6.38) 


Again, replacing a, b,c by ec, d?, e? respectively in (3.1) and taking y=ab, 
2—T, we obtain 


tr {a b (cd? e? + ce d?+@ece't+@ec+ed’c + e?c d?)} =0. (6.39) 


Further, replacing a, b, e by b, e, e? respectively in (3.1) and taking y=ad? 
and =I, we obtain 


tr{ad?(bee’?+ bec +cetb+cbhe?+ebec+e'cb)} =0. (6.40) 
With appropriate changes of notation, we see from (3.3) that 
trab@ee’?=0, trabe’cd#*=0, trad@bce?=0, (6.41) 
tra@bec=0, tra@ece’b=0, trad’cbe=O0. 
Using these relations in (6.39) and (6.40), we have, with Lemmas 1 and 2, 
trabed@e?+trebad@e?+ treab@e?+ trbac@e=0 — (6.42) 
ie trbcad@e’+ treba@e’=0. (6.43) 
From (6.42) and (6.38), we have 
tracb@e+ trbca@e’=0. (6.44) 
Interchanging d and e in (6.43) and employing Lemmas 1 and 2, we obtain 
tracb@e+ trabc@e=0. (6.45) 
Interchanging a and b in (6.43) and (6.45), we obtain 
tracb@e+ treab@e’?=0 


and (6.46) 
trbcad@e+trbacd@e?=0. 


From the relations (6.43) to (6.46), we see that each of the invariants formed 
from tr abed?e? by permutations of a, b, € is equivalent either to trabed*e? 
or to —trabecd?e?. It follows that, of these invariants, we may omit from the 
integrity basis all except tr abed?e?. It is easily seen from Lemmas 14 and 2 that 


trabced?=trebade?, 
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with analogous results for each permutation of 4, b, . Consequently, of the 
invariants formed from tr abed?e? and tr abce?d? by permutations of a, DB, €, 
we may omit from the integrity basis all except tr abcde’. 

Therefore, of the invariants formed from tr abed?e* by permutations of 
a, b,c, d,e, we may omit from the integrity basis all except 


trabcde?, trabdec?e?, tradcb’e?, trdbca'e’, 


trabee?d, traceb?d, trebead’, 
(6.47) 

tradebe?, trbdea’c?, 

tredea? b?. 


Collecting the results obtained above, we see that an integrity basis for the 
five symmetric 3 x3 matrices a, b, c,d,e under the full or proper orthogonal 
group is formed from the integrity bases for the five sets of four matrices which 
can be selected from a, b, ce, d, e, together with 

(i) trabcde, trabdec, trabecd, 
tracdbe, tracbed, tradbee; 


(i) trabcede*, treabde*, tredabe’, 
tracdbe?, trdacbe’?, tradbce? 


and invariants obtained from these by permuting a, b, €, d, e cyclically; 


(ii) trabecde’*, trbceade’*, ircaebde?’, 
trbcdeae*, trebaede*, treadebe?, tredbeae? 


and invariants obtained from these by permuting a, b, c,d, e cyclically; 


(iv) trabcd*e?, trabdc*e?, tradcbe?, trdbcae?, 
trabee’d’, traceb’d’, trebeadd?, 
tradeb’c?, trbdeadc?, 
tredea’b?. 
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1. Introduction 


In earlier papers, ADKINS & RIVLIN (1955) and ADKINS (1956a, 1956b) have 
considered the deformation of sheets and shells of a material which consists of 
incompressible isotropic elastic material reinforced with inextensible cords. 
Although the appropriate mathematical formalism can be developed for such 
problems, the presence of the elastic material makes it very difficult to solve 
any but the simplest problems. Furthermore, there are many situations of 
technological interest, in which the elastic material is either absent altogether 
or plays a secondary role in the determination of the deformational characteristics 
of the system. Accordingly, Rrvirn (1955) considered the problem of plane strain 
in a net formed by two families of straight parallel inextensible cords. 

In the present paper, we consider the deformation of a net, formed by two 
families of inextensible cords which, in the undeformed state, forms a surface 
not, in general, flat. The cords of the two families are assumed to form a curvi- 
linear net on the surface and the angles of intersection of the cords of the two 
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families are assumed to be independent of position on the surface. The radius 
of curvature of the surface formed by the net, in both the deformed and un- 
deformed states, is considered to be large compared with the intercepts formed 
by the adjacent cords of one family on a cord of the other family and this inter- 
cept is independent of the position on the surface and is the same for both families 
of cords. As in the earlier paper (RIVLIN 1955), it is considered that the cords 
are perfectly flexible, that the intersecting cords of the two families cannot 
move relative to each other at their points of intersection and that no two adjacent 
cords of a family are brought into contact as a result of the deformation. 


The deformation of such a net by forces applied to it is treated as a membrane 
problem, the membrane having two directions of inextensibility at each point. 
Appropriate governing equations for the determination of the deformation and 
cord tensions are derived in §§ 2 to 6. In § 7, we simplify these relations in the 
case when the net forms a right-circular cylinder on which the bisectors of the 
angles between the two families of cords form generators and lines of latitude. 
The cylindrically symmetric deformation of such a net is then considered in § 8, 
in the case when the applied forces are normal to the net in its deformed state 
and, measured per unit area in the deformed state, are constant over the net, 
and the edges of the net are displaced axially. In §9, we consider that the 
membrane is further surrounded by a smooth rigid coaxial cylinder, with which 
part of the net comes into contact when it is deformed. 


In §10, we consider a net formed by inextensible cords to be deformed in 
some known manner by specified forces and then subjected to a further infini- 
tesimal deformation. The governing equations for this infinitesimal deformation 
are derived. In § 11, they are specialised to the case when the initial deformation 
has cylindrical symmetry. 


In the problem considered in § 8, it is seen that the net remains undeformed 
under the action of the applied forces if the ends are held at a particular separa- 
tion. In §§12 to 14, we consider that the net is held in this manner and its 
edges are then subjected to small displacements. The manner in which the 


resulting small displacements throughout the membrane may be calculated is 
discussed. 


Finally, in § 15, we consider the case when small additional forces are applied 
to the membrane in the more general cylindrically symmetric state of deformation 
considered in § 8 and we discuss the manner in which the resulting small deforma- 
tion may be calculated throughout the membrane. 


2. The Equations of Equilibrium for a Membrane 
We shall take as our starting point the equations of equilibrium for a mem- 
brane in the form given by GREEN & ZERNA (1954). It is assumed that the 
equation of the surface in its deformed state is known, in parametric form, in 
a curvilinear coordinate system x, the two parameters being * X%, so that 


x =x! (X*) (2.4) 


x Throughout this paper we shall consider that Latin indices take the values Wi, 8 
and Greek indices take the values 1,2. 
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is the equation for the deformed surface. We shall denote by n%° the contra- 
variant stress components in the membrane in the surface coordinate system X 
defined by X!=const, X?=const. We shall assume also that the deformed 
membrane is acted on by tangential surface forces with contravariant com- 
ponents p* per unit area in the coordinate system X and by normal surface trac- 
tions p per unit area, the areas being measured in the deformed state. 

Let r be the position vector of a generic point P of the deformed membrane 


and let m denote the unit normal to the membrane at that point. Then, we 
shall define the quantities b,, by the equation 


dr. dn = —b,,d X*aX°. (2.2) 


a 


With this notation the equations of equilibrium are given (see, for example, 
GREEN & ZERNA 1954) by 


neh + ph 0 (2.3) 
and 
mn? bgt p =0, (2.4) 
where |, denotes covariant derivation with respect to X*, so that 
ap On™P | me a8 
Rabiner ON Be ate BI aes (2.5) 
where 
1 af 8p | O8ay a 
Li, =e -4 ee 26 
ty Be Oe DORE sax ue 


and g,, is the covariant metric tensor for the coordinate system X. The contra- 
variant metric tensor 2”? is defined by 


en" 2g = Of, (2.7) 


where 63 denotes the two-dimensional Kronecker delta. If y is an arbitrary 
three-dimensional rectangular Cartesian coordinate system, g,, is given by 


aye ay’ 
Sap axa axe en 
Defining g by 
& = £11820 — Sie» (2.9) 
we see from (2.7) that 
pile po./0; “2 = sie, and, gg" = — baie. (2.10) 


The components of the vector dr in the coordinate system y are given by 
dr =(22 axe (2.11) 
ax” 


The components of the unit normal n in the system y may be obtained in the 
following manner. We note that the two vectors with components (0'/0X1) dX? 
and (0y'/éX2)d X? respectively in the system y are both tangential to the mem- 
brane surface and are not parallel. Therefore, the vector with components 


ay! aye A 
fijk or gyn UXT AX? (2.12) 
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in the system y is parallel* to the vector mn. Since n has unit length it has com- 
ponents in the rectangular Cartesian coordinate system y given by 


ay? ay" ay™ os lite oyl za) (2.13) 


= lEpar Spmn axt axe Oxi OxX®| \"H* XT OX? 


Employing the relation 
25 ay Espn rm Orne AOD Orig (2.14) 
we see, with (2.8) and (2.9), that (2.13) may be re-written as 


1 ay! ay 
aan (ei aX! ax? } le 


/ 


We obtain, from (2.15), 


Cae ay! ayk B 
dn =(sx Ge Ej axl | aX"). (2.16) 


Introducing the relations (2.16) and (2.11) into (2.2), we obtain 


ar ay’ @ 1 dy? ay 5 
a8 ax ra “tik BI OX? |" a 


Equation (2.17) may be re-written as 


0 1 oy ay! ayF 14 dy dy? eryt 
bap = ijk { an a2 CRE . (2,48) 
PY OXP. | ph ORE OX OX gh @X1 ax? ax%axéb 
Since ay’ ayt ayr 
UP axe xt OMe 
equation (2.18) yields 
1 dyl dyk = ary? 
pe ee ne (2.19) 


gh aX! ax? ax*axe - 


Denoting the physical components of stress in the coordinate system X by 
tg, we have 


LENG o22\4 
m= (2 Nay, nt? = (EN tno, 
811 822 > 
12 git \t do at gas bY Cr 
a an | ica y o1- 
ta ~ tz a 
Introducing the relations (2.10) into (2.20), we obtain 
7a ae i ty te — (£11)! tos 
and 8811 &822 (2.21) 


gi pila gh en ote. 

The physical components /, of the tangential surface traction per unit area 
are given by 

h=(ei)tp! and fy = (G0)? p?. (2.22) 


Employing the relations (2.21) and (2.22) in equations (2.3) and (2.4) we can 
obtain the equations of equilibrium for the membrane in terms of physical 
components of stress and surface traction. In the next section this will be done 
in the particular case when the surface coordinate system X is an orthogonal one. 


* The positive directions for X1 and X2 are chosen so that (2.12) is in the direction 
of the outward drawn normal to the surface. 
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3. The Equations of Equilibrium when the Surface Coordinate System is 
Orthogonal 


If the coordinate system X is orthogonal, we have 
Sia — 0: (3.1) 
Then, from (2.9) and (2.10), we obtain 


Biren, gt = 1/g,, and g#?=41/g5. (3.2) 


o 
fe) 


From (2.6), we then have, with (3.2), 


r= 1 O81 = 1 og, eae AU fe Sa on Ogt, 
meer OX iz gh, C2” 2 2800 OX? ee: 
é 1 Og} 1 og 1 og3 
ee pee Shi 11 Les fe ore Sa 8 
21 281, OX* gh OX? 12, Ree Paes Soe OC gh, ox} G3) 
ee 1 0802 830 C83» and Bie eso, Chee 1 ogh, 
a 2g Oat 8%, ex} a i Dgee ON?) ph OX? 
From (2.3), (2.5) and (3.3), we obtain the two equations of equilibrium 
Ont ant? 2 Of, 1 est, ee 
axe axe ' Vor coxt "ph ox! 
11 22 
3 agt 1 083 gh, gh, 
2 SUL ee eee 2N dS 22 22 422 1 pl — 
(a “Ox? " gh, ox)" Cie, We 
and : (3.4) 
an* an®® CSL s Fe me ok tees ae 
axi ' ax? ga, OX? ge, Ox} 
22 ts! 
2a 1 esis) ee an 2 
; : 0 
aes ax T gh Oxt)” —~ g,, Ox?” Ge 
From (2.21) and (3.2) we obtain 
ti = O11, tee = 8a? and 42=ty = (211 S22)* 2. (3.5) 
With (2.22) and (3.5), equations (3.4) may be re-written as 
6(g2,¢ 6g? oge 
ee a ayia —0 
and Neer mst: Sag (3.6) 
py 2 dg3 1 
iatas f i Len +f aK tio + (811 S22)" fe = 0. 


Employing equations (3.5) and (2.19) to substitute for n** and by, in (2.4), 
we obtain, with (3.2), 


IR ES TACHA aS a 


t 
Ch EGE erp a (Bg ORO coe (611822)? OX1 OX? can Saag Gok i (3.7) 
+ (811 820)? = 


é 
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“a 


Equation (3.7) may be re-written as 


ay! ay! ax Oxf 
ETRE One ORY OX? 


ee 
x a ym 


ay 


r ax” ax” 


ar P (811820)? = 0 


\ 


E 


mL m 


hy 0? 4” 2he OP x ie toe ated 
ine BE year TL ERE | : F 
B11 0X1 0X (811822)? OX1 OX? Bae OX? OX? 


| 


Nh OM Oe 2h ee ee ae ee 


We ce Be ne TE Otay 
11 OXt OX (811822)? OX? OX? S22 OX? aX? 


and in this equation of equilibrium, as well as in the equations of equilibrium 

(3.6), Zag is given, from (2.8), by 

C. Le 
mn ava axp ? 


Sap = (3.9) 


where G,,,, is the covariant metric tensor of the coordinate system ~. 


mn 


4. The Deformation of a Membrane of Fabric— Kinematical Considerations 


We now consider the deformation of a membrane of the ideal fabric described 
in §4. We assume that in the undeformed state the distance between adjacent 
cords is # and that at each point of the membrane the cords of the two families 
are inclined at an angle 2x. We define the position of a point on the membrane 
in its undeformed state in the curvilinear orthogonal coordinate system X formed 
by the bisectors of the angles between the cords. The coordinates X% of a point 
of the undeformed membrane in this system are then the distances of the point 
from an origin measured along the coordinate lines. We choose the X? axis so 
that the cords of the two families are inclined at angles « and —« to it. 


Let us suppose that the membrane is deformed so that the point X* moves 
to x‘ in the curvilinear coordinate system x, the covariant metric tensor for which 
is G 

If dS is the length of the element in the undeformed membrane joining the 


points X, and X,+dX, and ds is the corresponding element of length in the 
deformed membrane, we have 


mn* 


GS\VeadX*d xX, = dX can 

and (4.1) 
(ds\*=e,,aX7ax". 

where g,g is given by (3.9). From (4.1) we have 


ds 2 5 
(as) ~ 80 L* 1! (4.2) 


where L* is the cosine of the angle between the linear element on the undeformed 
membrane and the direction of X* at the point considered. If the linear element 
coincides with one or other of the cord directions at X%, we have 


I’=cosa and L?=-+sing. (4.3) 
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Since the cords are inextensible, we have (ds/dS)?=1 for L* given by (4.3). 
Introducing this into (4.2), we obtain 


2 nz a 
mei Pig COS1K = Lon slnan — 1 (4.4) 
812 = 0. 


From (2.8), the second of equations (4.4) may be written as 


dy? ayt 

xt ext Ae) 
where y’ are the coordinates of a generic point of the membrane in an arbitrary 
rectangular Cartesian coordinate system y. Now, @y'/@X1 are proportional to 
the direction-cosines in the system y of the linear element in the deformed 
membrane which lies along the X? axis in the undeformed membrane and @y'/@.X2 
are proportional to the direction-cosines in the system y of the linear element in 
the deformed membrane which lies along the X? axis in the undeformed mem- 
brane. Thus, equation (4.5) and hence the second of equations (4.4) expresses 
the perpendicularity of these two elements. 

We can also define the position of a point on the undeformed membrane with 
respect to the curvilinear net formed by the cords of the two families. If 2% 
(a =1, 2) denote distances measured along the cords of the net which are inclined 
at angles « and —« respectively to the direction of X1, we have 


gi 2 @! 2) 
dn 5 (EE 4 P| and ae 5 (Se _ 
2 \cosa sin « ZeNCOSse sin % 
Whence, 
es «1 aoe ned eee. fat as 
ax! —2c0sa \eZ! ' 6Z? éx? ~~ 2sina \ éZ} oa 
and (4.6) 
ee COS & els SIN o - and pias eesg — sina — zs 
aZzt @xi | OX? eZ? ox? ax? * 
We may also re-write the second of equations (4.4) as 
ay" OV reg eee = 2 ii eitet ) 
eat (cosa SFr sin & Sy) = axt (COS % asa — Sina aoa). (4.7) 
With (4.6), (4.7) yields | 
oy? aye oF ay? oy (4 8) 


@X1 OZ! ~~ AX" OZ? 


We now consider three elements of length at a point of the undeformed membrane. 
One of them lies along the X! direction and the remaining two have the same 
lengths and lie along the Z1 and Z? directions. Equation (4.8) expresses the fact 
that in the deformed membrane the first of these bisects the angle between the 
remaining two. 

We thus conclude that the curvilinear net on the deformed membrane defined 
by X1=const, X?=const is an orthogonal net and at each point of the deformed 
membrane the coordinate lines bisect the angles between the cords. The directions 
of the coordinate lines are therefore principal directions for the deformation. 
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5. The Physical Components of Stress in the Deformed Membrane of Fabric 

If t, and t, are the tensions in the cords which were initially inclined at 
angles « and —« respectively to the X? direction, we can easily show (see, for 
example, ADKINS & RIVLIN (1955) or RIVLIN (195 5)) that the physical components 
of stress t,g are given by 


g2, cos? % g2, sin? a 
A= Ee (t+), Te2—=—™ ae (tT + Ta) 
839 it (5.1) 
__ sin «cos 
he = toy iro (t — Te) 


6. The Equations of Equilibrium for the Deformed Membrane of Fabric 


The equations of equilibrium for the deformed membrane are given by (3.6) 
and (3.8). Introducing the expressions (5.1) for the physical components of 
stress into (3.6), we obtain 


ser Leia 15) COS" Oo on [34 (t% = T.) | sin «cosa + 


0g}, eh, 083, 2 h 3 
SP as ax? Gt 62) SUE POs OxT (v1 + Tz) sin? + A (841 820)? ff =O 
1 


and (6.1) 


é a 
fee [g3e(4— T) | sin 4 COS % Si EEO [g30(t% an T.) | sin? - 


gh, : et, og}, 2 3 
ax (t — T.) Sina cosa — + Be (t, + t,) cos?a + h(g,4 800)? fo =O. 
522 


These equations may be considerably simplified by making use of the first of 
the relations (4.4). Differentiating this relation with respect to X1, we obtain 


1 Og2, 3 
et, a cos? a = — gh, —82? sina. (6.2) 
Introducing this result into the first of equations (6.1), we obtain 
0 4 
“OxXL [213 (t%1 + T,) | cos? + —— a [e13 (tT — T)| sin a cosa + hg, gf, =0. (6.3) 


Again, differentiating the first of equations (4.4) with respect to X2 and employing 
the result so obtained in the second of equations (6.1), we obtain 


4) 0 r) x 
Oxi [go0(ty ra T.) | sin % COS % ++ Sa) eevee -|- T.) | sin? o + h 85083, (gO. PU(G4) 


Employing the relations (4.6) in equations (6.3) and (6.4) we obtain 


a (a) 
yar (811 %1) + oz? (G11 T) + G11 30/1 =0 
and (6.5) 


a) a) 
Oz (a2 T%) — gze (82272) + <n baabts fs = 0. 
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Introducing the expressions (5.1) for the physical components of stress into 
(3.8), we obtain 
where 
ay! ay* ax? ax4 


P= Es; = = a 
i} oa? ant OR ox? 


| Geigy 2 Olan meee CP 
a OOM HE DS ~sing cosa + in? «| 4 
een aX) QX1 Tees cei eee @X®exe | 
ey? On” Ox” Cie OKs Cn” On 
cos? ++ 2 —_— ——_ sina cosa'4 -— sin? 
" oa a4” | axl axl He 2 srata yea Se COSH thas a Clot alt, 
Oe abyhnaah ext 08 a 
18 oy Ox4 OX! aX? 
| CL 5 i ar pees hed ie 
aoe — COS? a” — 2 — —~ sin % COS 2 |e 
er @X1 0X! QXi 6x2 1% COST Fea oe IM a 3 
eryt Ox” ax sg! Bree Bo! Oude 
| COS) = 2 pe ES ITNOACO sealer 2 
' 0x" 0x" | OX! aX} pakke as axe 6x2 1 
and 
R=hhP &1 809. 
We have 
ay? dy! ay i . 
Sik Gem ax? axl MPO] axl oe 


Employing this relation and the relations (4.7), we may re-write (6.7) as 


P Oy | Oxmmo dt 084 
= Enpq al 
Ox OX Nem OLL OZ 
Oyi ay” ax? axl ayt ax” ax” 
see thes pies p n ? 
CERO aOLX Ss OKA MO GON ay OZ ARO Z- 
6.9 
Gas Oy" | Ox? 0x2 029" 6) 
SPT lax! | 0X8 aX OZ? 672 * 
tag ay! ay" OxP axl ety = ax™ ax" 


118 aye Ox? OX1 AX® ax" ax" AZ? OZ? 
and 
R=hP 11 820- 


7. The Governing Equations for a Cylindrically Symmetric Deformation 


We shall now assume that in the undeformed state, the membrane forms a 
right-circular cylinder and that the lines X?=const are circumferential lines on 
the cylinder and the lines X!=const are generators. We assume also that X?=0 
is the equatorial circle on the membrane. We shall take the coordinate system 
x’ to be a cylindrical polar system 7, #, z with the z-axis coinciding with the axis 
of the membrane, the origin lying in the plane of the equatorial circle and ? 
measured from the axial plane through X1=0. 
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The covariant metric tensor G,,,, for the coordinate system is then given by 


Ose © 
|| Gen | = || 0, re, Oj]. (7.1) 
OF. Cen 


We now assume that in the deformation a particle at X* on the undeformed 
membrane moves to 7, #, z on the deformed membrane, where 


y= t(X*), Fe Ma “and “2=—2(X?): (7.2) 
Then, from (3.9) and (4.1), employing the notation X?=2’, we have 
fii = (ra)? and g5_ = (dr/dX’)*+ (dz/dX)? = (do/d 2)’, (7.3) 


where o denotes the curvilinear distance of the point considered on the deformed 
membrane from the equatorial circle, measured along a generator. 
From (4.4) and (7.3), we obtain 


do 


nan =f * sinta =“ (7.4) 
a (7x) sinta =4, 


and we note that the second of equations (4.4) is automatically satisfied. 
We shall also assume that the applied force system is such that f,=0 and f, 
and # are functions of 2’ only. Then, from symmetry considerations we have 
py Sh ay), (7.5) 
which, with (5.1) and (7.3), yields 


ieee 2T 7 COS? % i 2ta(do/d2) sin? « 
i) Sap daase l=" hy 


and -2,.=0. (7.6) 


Introducing equations (7.3) and (7.6) into equations (3.6) and bearing in mind 
that /,=0 and, from symmetry considerations, t is a function of X? only, we 
see that the first equation of equilibrium is automatically satisfied and the second 
takes the form d his 
“ga has) => — tage? fg =O (7.7) 


In the particular case when /,=0, 7.e. the forces applied to the membrane are 
normal to it in its deformed state, we can re-write equation (7.7) as 


d 
oie (7tg9) —hy =O. (7.8) 


The remaining equation of equilibrium may be obtained from equation (3.8) 
by introducing the relations (7.2), (7.3) and the last of the relations (7-6): In 
order to simplify the analysis we choose the rectangular Cartesian coordinate 
system y in such a way that the y!, y?, y3 axes have radial, azimuthal and longi- 


tudinal directions respectively at the point of the deformed membrane considered. 
Then, 


Membrane of Inextensible Cords 457 


Also, 


and, otherwise, gh yi (7.10) 
OPEC HZ 


In this manner, equation (3.8) becomes 


A ee tog e d*z dz aa) =? 


ne | Es su! . = = 
y do %1' (dojdd)?\do dX? do ad? (7-11) 


Since the deformed membrane has cylindrical symmetry, the lines of latitude 
and longitude on the membrane are principal directions of curvature. Denoting 
the curvatures in the latitudinal and longitudinal directions by x, and x, re- 


spectively, we have 
dz 


1 Gel az 
ii ee and 45:s3 ae (35 : (7.12) 
Introducing the relations (7.12) into (7.11), we obtain 

Ah 1 + Hala p. (7.13) 


Equations (7.8) and (7.13) are, of course, the well-known equations for the 
equilibrium of a membrane having cylindrical symmetry and deformed in a 
cylindrically symmetric manner by a cylindrically symmetric system of forces 
directed normally to it. 

By substitution for #,, and ¢,, from (7.6) in (7.8) and (7.11), and employing 
the relation (7.4), we can re-write equations (7.8) and (7.11) in the forms 
aia tT = const 

d (az) me (71.4) 
an 

(P+Q)t+R=0, 


where 


PA Ge GY We WEN se. (Ga ON ea 

lila fea ae ds? d= as? LS pe 

and ane do\? (7.15) 
ed i ) ; 


These equations can also be obtained directly from (6.5), (6.6) and (6.7) by 
employing the relations (7.2), (7.3), (7.5), together with /,=/,=0. 


8. A Simple Case of Cylindrically Symmetric Deformation 


In this section, we consider that the membrane initially forms a cylindrical 
surface of circular cross-section and is maintained in a cylindrically symmetric 
deformed state by a uniform normal force # per unit area, measured in the 
deformed state, its edges, which initially lie in the planes 2 = +/, being held 
in the planes z= +c and with radii 7,, as shown in axial cross-section in Fig. 1. 
We denote by 7, the radius of the equatorial circle in the deformed state and 
by 7 the radius of the circle formed by the points on the deformed membrane 
for which dz/dr =0. 
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From equations (7.6), we have 


tga , Oe Vas 84 
ban yf (a5) eel es ( ) 


Employing (8.1) to substitute for f,, in (7.8) and (7.13), we obtain 


d a? (do \2 9 
: Be =f, 
ay (0-7 (az) x ss 
and (8.2) 
GAG Nea ee 
ba +m a ae) tan? a} 2, =P. 
42 
| 
[t----- 
| 
| 
rs pa 
| 
| 
| 
| 
Fig. 1 


Eliminating ¢,, from equations (8.2), we obtain with dp/dX'=0, 


r( 2) sin? COS” % 
d pe Ne ol eee (8.3) 
dy | r ie ag. GO \? ae vi y? - Pa Sn ‘ 
Hg COS a +, ( G5] sin? « ae costa + (45) sin? « 
With (7.4), and employing the notation 
r(1— 5 cos? «| 
p=", + (8.4) 
#X_+ (x, — %2) = COS? w 
equation (8.3) may be re-written as 
SS 2 
io 3 COS? & é 
dy v2 E < (8.5) 
1— —, cos? « 
Equation (8.5) may be solved for @ to give 
A 
YP Tro Tune) ) (8.6) 


ye a 
(1 —-, cos? a] 
a 


where A is an integration constant. From (8.4) and (8.6), we obtain, with the 
expressions (7.12) for x, and x5, 


Aa a Pe eee cd y 2. 22 1 re 2,\% 
[ga 8 a agg) gk ON ge gt ee ena ait 
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Solving this equation for dz/do, we obtain 


az 4 & 3 
=(s4 v4 B) (1 a ae cos a) ; (8.8) 


where B is a constant of integration. Introducing into (8.8) the conditions that 


De eee (ye when if = 14) 
dy do 
and de (8.9) 
7 a 1 when 7 =, 
we obtain 
( 54 Ya B) (1 us cos? a)! = 
and (8.10) 
1 2 B\(4 vr 2 bp 
(4 1; )( Ge 008 a| = 0. 
Solving equations (8.10) for 2A and B, we obtain 
2 2 
B=— 4 = _ (8.11) 
oe (72-7?) pees: cos? a)" 
From (8.8) and (8.11), we obtain 
dy 1 2 a Me 
ce alee eg pale 752 (1 - cos? a)| (8.12) 


In this and the following equations, the upper sign is applicable to the region for 
which 2 is positive and the lower sign to that for which z is negative. From (8.12) 


we obtain 
rf 


eh . 81 
fi nay (1% cost a) (8.13) 


Yo 


the lower limit of the integration being taken as 7 in order to satisfy the condition 
g—0, when 7 =75. 

The equation, in the cylindrical polar coordinate system, describing the form 
of the deformed membrane is given, from (8.12), as 


x 1 


ae (v2—77) (1 us cos? a) dy 
Vike *{ ad (8.14) 


2 4? 
1 we (F 7)? (4 5 cos? a)| 


Yo 


the lower limit in the integration again being taken as 7, in order to satisfy the 
condition z=0, whenr=7. 
If, as we initially assumed, z= +c when r=7,, we have, from (8.14), | 


ts 4 


(v2—73) (1 és cos? a) dy 

2A t a 

— / i fs (8.15) 
‘ 4 —sge ri)? (1 cost a)| 
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From equations (7.4) and (8.12), we have 


2 2 4 y? 3 
oS ==> — (1 — os cos? a) Retr) hen ie (1 — =, cos? a. , (8.16) 
whence 

; d 

=> sin« | — a aoa eg s 3 ars, (8.17) 

(1 me cos? a 1 4A (v2—r3)? (4 = cos? x)| 

Since / = +1, when 7 =7,, we have from (8.17) 
; ay 

t=sina | A ; ; - x Sere (8.18) 

(t= qe CO” a| 1 7H (v2 v2)? (1 mT cos? x)| 


i) 


If 7,7, and 7, are specified, then equation (8.14) can be used to determine 
the shape of the deformed membrane and equations (8.18) and (8.15) may be 
used to determine the length of the undeformed membrane and the positions 
in which the edges of the deformed membrane should be held. 

If, on the other hand, c,/ and 7, are specified, then equations (8.15) and 
(8.18) provide two equations for the determination of 7) and 7,. Equation (8.14) 
may then be used, as before, to determine the shape of the deformed membrane. 

If 7) and 7, are known, then from (8.8), with (8.11), we can determine dz/do 
at each point of the deformed membrane and, from (7.12), the principal curvatures 
can be determined as 

1 


eS y 
: Wee wy (7? — 73) (1 qe COS & 
an 
mt =35 (1 — © cos? a) * es (372 — 7?) sos 
ae A a® i a ui" 


From (8.2), (8.4), (8.6), (8.1) and (7.4), we obtain 


Y 
A pcos? a — 
Ai Uh mr) g 
(1-7, cos? a| 
and i (8.20) 
ee: , 
as yr By 
r(1—7 cos? a} 
a 


From (7.6), we obtain, with (8.20) 


h Abh 
a eae: (411 t2)* = oie 


(8.21) 


3 y? ’ 
2asin« (1 ——~ cos? a| 
a2 


HOFFERBERTH (1956) has derived a formula analogous to (8.14) in the case 
when the cords are extensible. However, his expression for z takes the form of 
an integral in which the integrand contains the cord extension whose dependence 
on 7 is unknown. In the particular case when the cords are inextensible, the 
extension is zero and HOFFERBERTH’S formula reduces to (8.14). 
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9. Further Development of the Problem 


We now consider that the cylindrical membrane discussed in § 8 is located 
coaxially with a rigid circular cylindrical ring of radius 74, so that when the 
membrane is deformed by the uniform normal force #, its edges being held at 
¥=1,,2=+¢, it still assumes a form having cylindrical symmetry, but is partly 
in contact with the ring, as shown in axial cross-section in Fig. 2. We shall assume 
that there is no friction between the membrane and the restraining ring. 


ed 

| | 
Bonaplagepe™ 
| 

oat ie a 
| 

| / 

r 70 

| 

| 

| 

| 
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Fig. 2 


For the region MN in contact with the restraining ring, we have 
jee Aly, And) “x, ==. (9.1) 


From (7.4) we then obtain, on MN, 


From (9.2) and (7.6), we see that on MN 


9 COS? w sin 2ta %o2 ee 
em 2T 1% ad =. AINGs be fon = (1 = cos* a.) sin a. (9.3) 


12 4 hy, 
ah (1 — —_ cos? a) 2 
a 


On the portions QM and RN of the membrane which are not in contact 
with the restraining cylinder, equations (8.1), (8.2) and (8.3) still apply, so we 
obtain, as in § 8, the result (8.8). It is easily seen, from NEwron’s third law, 
that at M and N the membrane is tangential to the restraining cylinder. We 
therefore have 

dzidy =O when r=”, (say) 


and 3 (9.4) 
adzidg =1 when r=7%.- 


These relations replace (8.9) as the appropriate boundary conditions for the 
determination of A and B in (8.8), and we obtain 


B ihe Ne a (0.5) 
2A \ 752 = k 
(79?—72) (1 Sih cos? | 
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If the z-coordinates for M and N are +2, we obtain from (8.8), together with 
(9.5) and the fact that dz/do=1 on MN, 


ee dy 
; sar f RAE a 
I 7A2 (v2?—77) ( me oS a)| 


and (9.6) 
M 1 Ve 2 
Ey (7277) (1— Z cos? | dy 
A ee ee 
1 2 y2)2 (1 é 2 »\ | 
4 4AB (72—73) we 088 x) 


the upper signs being applicable in the region QM and the lower signs in the 
region kN. 


From (9.2) and dz/do =1 on MN, we see that on MN 


Thus, 
PIS a (1 = fo cost a) z% at M. (9.8) 


On QM and RN, equation (8.16) is still valid, so that, with (9.8), 


: 7%" se, 
E—sin a(t — 2, cos*a) 2 
az 


(9.9) 
= : dy 
: (1-5 cos a| |1— 4A2 (Gon) (ta cos a)| 
We obtain immediately 
12 —_ 
1 = sin «(4 — “4,- cos? ea = 
a2 
: (9.10) 
sina f _——_, wiay ——--—, 7 
Y 2 1 . y* i 
: (1-7 cos? x| 1 4A2 Ge #)? (4 a cos? a) 
Also, from the second of equations (9.6), we have 
" 1 2 4 
oer (v?—7?) (1—— cos? a) dr 
C2 = usar 
0 7 7 tr: (9.11) 
f 4A? (v2—4})? (1 = cos? «)] 


Equations (9.10) and (9.11) provide two equations for the determination of 
7 and 20 if c, 1,7 and 7, are specified. Equations (8.19), (8.20) and (8.21) are 
still valid for the determination of #1, fy9, t, %, and x, on QM and RN. 


Membrane of Inextensible Cords 463 


10. The Superposition of a Small Deformation on a Deformation of the 
Membrane 

We now consider that the membrane discussed in § 6 is first subjected to 
a known deformation in which a generic particle initially at X% in the coordinate 
system X moves to x’ in the curvilinear coordinate system x and surface tractions 
with physical components /, and #, measured per unit area in this deformed state, 
act on the membrane. /, are, of course, the components in the directions of 
the parametric curves X?=const and X!=const on the deformed surface and 
p is the component normal to the deformed membrane. The membrane may 
also be acted upon by edge tractions or subjected to edge displacements, but 
we shall not specify these in detail. This state of deformation will be called the 
first deformed state. 

We consider further that the membrane is acted on by additional small forces, 
so that the generic particle undergoes an additional small displacement ew’ in 
the coordinate system x. The resulting state of deformation will be called the 
second deformed state. We shall employ the notation /* and p* for the physical 
components of the surface traction in the directions of the parametric curves 
on the membrane in the second deformed state and normal to it. 

These are, of course, measured per unit area of the membrane in the second 
deformed state. We shall write 


fe =fot+ef, and pt=p+ep. (10.1) 


We shall assume throughout our analysis that ¢ is sufficiently small so that 
terms of first degree in ¢ may be neglected in comparison with terms of zero 
degree, terms of second degree in comparison with terms of first degree and so on. 
In this section we shall derive the differential equations determining w'. 

Writing 
A Ee, Gry — Cy, 1 Gog (10.2) 


where G*,, is the covariant metric tensor for the coordinate system x at the 
point x;*, we apply equation (3.9) to the second deformed state of the membrane 
by replacing x’ and G,,,, by x*' and G*,, respectively. Denoting by g*, the 
expression for g,, then given by (3.9), and writing 


Bre = Sup te oan. (10.3) 


we obtain, by equating coefficients of ¢°, the expression (3.9) for g,, and, by 
equating coefficients of ¢, the expression 


ap a Gaal 


(10.4) 


GAn OU, OR OWN On" ORF G 
CxO CR OK OX" J. | OX" OKe 
Applying equations (4.4) to the second deformed state of the body by replacing 
&ag by gx, and proceeding in a similar manner, we obtain equations (4.4) and 
811 C0s*% + 8osin*x=0, g12=0. (10.5) 


Equations (4.4) express the fact that the cords are not extended in the first 
deformation and equations (10.5) express the fact that they are not extended 
in the second (small) deformation. 
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We now consider the equations of equilibrium (6.5) and (6.6). Employing 
the notation 1,, t, and t*, r¥ for the tensions in the cords of the two families 
in the first and second deformed states respectively, we write 


c*=t%t+et, and t =1)+ &7p. (10.6) 


Then applying equations (6.5) to the second deformed state of the membrane 
by replacing £11, 22, TM, T2,f, and fy by gti, a2, Ti, Ta, fr and fz respectively 
and making the substitutions (10.1), (10.3) and (10.6), we obtain equations (6.5) 
and the two equations 


= 0 ~ = é a 
a ($1171) + 372 (S11 T2) + EA (T1841) + OZ? (T2811) + 


h 811 1 S22 | h == (i) 
T cosa ihe e i 2 Boo | ‘) (10.7) 
and 
@ - é = é g z 
OZ) (8227) — Be (S22T2) + zi (T1 G22) — OZ? (T2822) + 


0 Garey 4 Brat 4 ht) = 9. 
sng Oth tl 2 et geass 


sin « 


Again, in (6.6) we replace P,Q, R, %], tT, by P*, Q*, R*, tf and tz for the 
second deformed state, where P*, Q* and R* are defined by the expressions (6.7) 
in which x’, y', p and g,, are replaced by x*', y*', p* and g%,, where y*' denotes 
the coordinates in the rectangular Cartesian coordinate system y of a point at 
x*' in the curvilinear coordinate system x. Then, with the notation (10.1), 
(10.2), (10.3), (10.6) and 


Pt=PteP > OF 20+ 20, Ka koe (10.8) 


we obtain equations (6.6) and 


PHO + YP wBO RR = 0: (10.9) 


11. Small Deformation Superposed on a Cylindrically Symmetric Deformation 


We shall now consider that the membrane has the form of of a right-circular 
cylinder of radius a in its undeformed state and that the bisectors of the angles 
between the cords lie in the circumferential and longitudinal directions on the 
membrane, as in § 7. We shall also assume that the first deformation is cylindri- 
cally symmetric and is described by equations (7.2) in the cylindrical polar 
coordinate system 7, ,z, and is maintained by a uniform surface traction p 
per unit area applied normally to the surface of the membrane in its first 
deformed state. 

The fact that the cords are not extended in the first deformation is then 


expressed by equation (7.4) which, with the notation dr/dX =r' and dz/dX =z’, 
may be re-written as 


2 
Costa, + (7/2 + 2°2) Sin7 =="). (11.1) 
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The equations of equilibrium for the first deformed state are given by (7.14), 
which may be written as 


ig 22) = const (£422) 
and 
2P1+R=0, (11.3) 
where 
Nees - (3° 7’ — 7" 2") sin? o — 9 2OCOs O 
and a (11.4) 
R Ship ey ae 2) 


and 7 is the cord tension in the first deformed state. 
Taking x'=r, #, z and employing the notation u'=vw, v, w, we see from (10.2) 
and (7.1) that 


Gee Ory and G0" (nn-= 2 2). (11.5) 


From (10.4), (7.1), (7.2) and (11.5), we obtain, with the notation X!=@ and 
DG 


Se 2 ov, u = ( , OU , OW 
= a ===) ae Estie oe 
§11 A (r 260 | ak £22 es 2 al 


and (11.6) 


= chap Eig ov , OW 
Siz 00 a 0a 


Introducing (11.6) into (10.5), we obtain 


nib (# +735) cos?a 4 [7 hed f3! a) Sin?et oh 


q a\a XC) ha 
a RO hte CON, OW) & a7) 


0 ' a 6: '* 30 


Since the first deformed state is maintained by a uniform pressure p per 
unit area applied normally to the membrane, we have 


ae 
and vidernh (11.8) 
Ty — T9 =T. 
With (11.8), (7.3), (11.6) and (4.6) and, for convenience, using the symbols 2, 
and Z, in place of Z1 and Z?, equations (10.7) become 


2 OG OT QU = Beat 
Blan, ! az) ees aay oe) eet 
at 
Cee | Walla me CO. RY) ns Le a 
eles beng 56) cose + ar a a) ae ae 
d 9 
= ia" 2) (3 a t 2(7' 7" + 2° 2"") (% +7.) sina + 
OL; OZ. 
7) , OU 7 OWN .- , OU pen Bas 
4t | as 2 oa) sina + 4(r aa t# Set sin a ++ 
+ h(r'? + 2’) 4 rus 0, 


where the notation t’= dt/d& is used. 
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Again, we obtain from (10.9), with (11.8) and (41.4), 


t(P+0)+ “\(2 "pm FB) sin?a — 7, 2' costa! 2a) =0, (11.10) 


where, from (6.7) and the definition of P,Q and R given in § 10, we obtain 


= QVz O2u Pa Cru 
P+Q= a (co 20 BOF + sin ass 
avr ew rw 2 2Y > aaetN 
ee, [costae sin2 % + z (r" sin? x ee cos a.) + 
+2 oe z ar sin? « — * cost a al. 
19 INE De BD US ae ) (= a i ee Sp ay) 2” sin? a 
eb a PE: Sie ae ao) 
and (44445 
ota y2 , Ou , Ow | if 12 12 iu BCL ath 
R=2hp|7(r ae tt oe) tot +2'9 (= +7 4 
thp es (7/2 + 22) 
a> 


12. Small Deformation of a Right-Circular Cylindrical Membrane 


Let us suppose that in the analysis of § 11, the first deformed state coincides 
with the undeformed state, but a uniform normal surface traction # is applied 
to the membrane. It is easy to see that this is a possible state of deformation 
if the edge conditions are appropriately chosen. Equations (7.2) become 


t=4, (Ph =Oe and (12.1) 


so that equation (7.4) is satisfied. The first of the equations of equilibrium 
(7.14) is satisfied if t—const and the second is satisfied if the value of this 
constant is so chosen that 

ahp 


ag 2cos2« ° (1 


2.2) 
It is thus apparent that a possible state of deformation of the cylindrical membrane 
is given by (12.1) if the edge displacements are zero, even if +0, and that the 
tensions in the cords are then given by (42.2). From (7.6) it then follows that 
the edge tractions must be longitudinal and given by ap tan?«. 

With (12.1), we obtain from (11.7) 


02) 


) (12.3) 
ssn ae 0s 


(= +a 
and 


Again, from (11.9) we obtain, with (12.1) and (12.2), 


OU ai 1 hi — 
ae OZ, 4r(2 oe oa Ty gr) Cos + oe oF —()) 
(12.4) 
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Also, from (11.10) and (11.114), we obtain 


t(P +0) == (Z, +7,) costa +R =0, (12.5) 
with 
Mee Ty U ty OU Ov 1 dw 
P+ Q =2(costa 2a, + sinta 2) 2(2 ay os — pyr) COS? 
and (1 2.6 
ss u ov Ow 
B 2np (2 ! XC) ! = e 


13. Small Deformation of a Cylindrical Membrane with Specified Edge 
Displacements— Governing Equations 
We now consider that the second deformation results from specified dis- 
placements of the edges of the shell at 2 = +/, and that f,==0. Equations 
(12.4) then become 


Oly § Obs | ( ou |, a?” = 
and - = (13.4) 
OT, OT, | Tee ay 0 
Wa te od? ; 


and, in (12.5), P+ Q is given by the first of equations (12.6) and R is given by 


D U du, ew 
R= 2hp(% +a 20 +35). (13.2) 


We can eliminate uw, w, tT, and tT, from the five equations (12.3), (13.1) and 
(12.5) in the following manner. 


Differentiating the first of equations (12.3) with respect to 0 and employing 
the second of equations (12.3) to substitute for 0w/00, we obtain 


a Ce GOs 
ae 2 2/_ 2 2 
COS” & a | 32 COS" a ye sin a) 6 (1 3.3) 


Employing the relation (13.3) in the first of equations (13.1), we obtain 


b Ata =0. (13.4) 


OT, ea ev sin? « dt zw Eee 
potas aX? cosa ax? 
4133- 
ie Oey Dea ev sin? a ey rw ae Ce) 
6Z, 0X? cosa ox? j 


Differentiating the first of equations (13.5) with respect to Z, and the second 
with respect to Z, and adding the equations so obtained, we derive with (4.6), 


o2 
Lda 


= 
as \ “gees " om3 


= s ; 2 Ort 
(t+ To) =47 sin? a, gee 7) 


(13.6) 
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Now, differentiating (12.5) with respect to 0, we obtain, with (13.2), 


P t a oe C20 Ve 
(35 Bee cota “(SB 4+S 3) +2 hela oe. aor 3025) es, 


Substituting in (13.7) for @u/@O from (13.3) and for @w/0@ from the second of 
equations (12.3), we obtain 


= 


a) 1 OT. a 2hPpa |... ayy Ce 
ee Ne * costa (52 Fai ae (sin? « — cos? a) 352 = 01 (43.0) 
Now, from eu 
= i a { au ‘ CU es 
- (P+Q) = 245 (5 FICE COS?H + Axe sin a) — on 
4(2 ou 4 3 oy 1 ew ) costa. j 
& 00 ~~ 00? ' a 0002 


Substituting in (13.9) for du/dO from (13.3) and for ¢w/é@ from the second of 
equations (12.3), we obtain 


el eee 
ao EY) 
ot EO eae ey ee 
ate (a8 cos? — ag Sin | aor OS &  aeqisin a) v (13.10) 
— 2 costa ig, oak (cos?a — 2 sin*a) . 
00? od? 


Forming the second derivative 6?/0Z,0Z, of (13.8), introducing the results 
(13.10), (13.6), (42.3) and (12.2) and bearing in mind, from (4.6), that 


Oe oct egy Gee (43-11) 
Gz ant sor as?’ 
we obtain 
Mie We a price ay e 
la (cos Cyc sin a aoe) (cos* a a + sin? « 05 | - 
(13.12) 


+ costa (cos? « x Sh cinte oo ~~ eee pee Oe 
00? 0} \ 002 OX 


14. Small Deformation of a Cylindrical Membrane with Specified Edge 
Displacements— Solution of the Problem 
If we exclude the possibility of dislocations or tears in the membrane, wu, v 
and w must be continuous single-valued functions of position on the membrane. 
They must therefore be continuous periodic functions of O with period 22a. 
It follows that we may express uw, v and w in the forms 


co foe) 
u => Wee Vv = 2 ven 
n=0 


and t tip (14.1) 
w= 2 Wy, oe 
n=0 


where w,,, v, and w,, are complex functions of & only. 
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We therefore consider a solution of (13.12) of the form 
jis ye” Oi: (14.2) 


Substituting from (14.2) in (13.12), we obtain 


& tan a nyt) + (= tan a) (n? — 3 cot?a) vl”) + 
a e if (14.3) 
eS (“ tana. G — cot?a — n*) Uy ate (1 =n?) Uv, = 0, 


where v', v”) and vj’ denote the sixth, fourth and second derivatives respectively 
of v,, with respect to 2. 

Corresponding to any solution for v of (13.12) of the form (14.2), we can obtain 
equations for w and w from (12.3). Substituting from (14.2) in (12.3), we see that 


in Ola 


U =U, € AN oe =e ee (14.4) 


where uw, and w,, are functions of ' only, satisfying 


( aE ine) cos?a + w, sin?« = 0 
and (14.5) 
a Vp nel? ae Lie =i), 
These equations yield 
2 n az 
u, = —ina (*, — tan? « v, and wy, =7 = ee (14.6) 
unless =O. 


The expressions for 7, and 7, corresponding to a solution for u,v and w of 
the form given in (14.2) and (14.4) may be obtained by taking 


Z, =T eile and =F, = 7) en 1, (14.7) 


n (43.5) and (12.5), where 7) and 7?) are functions of S’ only, and employing 
the relations (4.6). From (13.5) we obtain 


in sin? ~ 
—* cosa t® + sina t” = — 2rav, — 2Tw,, sina 
a COS &% 


and (14.8) 


an = : —(2)’ rr Sin? 
—— cos at?) — sina tT?) = — 2tav 7 + AT Wn sin @. 
a 


From (12.5), we obtain, with (12.2), (12.6), (13.2), (14.2), (14.4), (14.7) and (14.6), 


1 costa (7) + # ”) =hp(— “tt, +atanta uy — inv, + 0] 
a 
a on at i (14.9) 
=inhp le Lee Un tan’ « ee 
From (14.8) we have, with (14.6) and (42.2), 
a Tipe. 21ehpsina ow 
=* cos a (7f)— tp )+- sin o (79 + 72”) = —47w, sina= SE (14.10) 


470 R. S. RIvyin: 


Substituting for 7!) + 7 from (14.9) in (14.10), we obtain 


ver 


a* a* 
(n?—1) v, +3 ae ays tan‘ « v!”) 


_ (14.44) 


1 costa (7 — 7) = —ahptana 
a 


Equations (14.9) and (14.11) may now be solved for ti) and 7,”) yielding 


2) me 7 a2 1 a* 4 (iv) 
a costa tT =inhp |(n?— 1) Vat ag Oye ge FOO, ae 


2 i a ad a t 4 (v) 
—ahptanae \(n —A1)v,+3 a Un — ye tanha vy 
and (14.12) 
2 =(2 : a2 + a* 4 (iv) 
— COS? 4 Ty” =inhp | (n? 1) ¥, +5 Un — —@ tanta v, | + 
/ a? wr a 4 (v) 
+ ahptancs| (n? 1) 0, +3 ety ee tan*« v,, . 


We shall discuss here the solutions of the problem for m=O and »=1 and 
leave discussion of the more general case for a later paper. 


Case n =0. Equation (14.3) becomes 
aul) — 3 cot4a vi) = 0. (14.13) 
Also, the second of equations (14.5) takes the form 
Yo =0, whence v,=const, (14.14) 


which automatically satisfies (14.13). The remaining equation (14.5) now takes 
the form 


“* cos*a + wy sin?a = 0. (14.15) 


Also equations (14.8) take the forms 


TY =—27tw, and 72’ =— arm. (14.16) 


From (14.16) we obtain 
tT =—2rm+A, and 7 =— 27H) +4,, (14.17) 


where A, and A, are constants of integration. Also, introducing the result (14.14) 


and the expressions (14.2) and (14.4), with —=0, into (13.2) and the first of 
equations (12.6), we obtain 


/ 

D 7a) : U W, 
P+Q=2sin?au) — 2 (2 0 ale =) cos? % 
and a ie 


R= 2hp (%® + 1%). (14.18) 


Introducing (14.18), (14.17) and (12.2) into (12.5), we obtain 


+ 9 ” 3 2 . eeeA 
Sino Wee gq COS & Wy = cos? a, (14.19) 
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where A =A,+A,. From (44.15) and (14.19), we obtain 


: wt A 
sint a ug — = Cost a Uy = 
a 2 


cos? % sin? «. (14.20) 
at 


The general solution of this equation is 


Ue A pete) A. Degen) Zia oe tan? «. (14.21) 
Introducing (14.21) into (14.15), we have 
MW = — aa Cola Agel A122 Be (Herero aie) oe o : (14.22) 
Integrating equation (14.22), we obtain 
= oA A, el¥3eot*a) Fa aa Bye7 (13eot*#) Sa 4 a Za, = (44.23) 


where & is a constant of integration. 


As an example, we consider the specific problem in which, in the deformation, 
each point of the edges X’= +/ is displaced radially by a constant amount ¢U 
and parallel to the cylindrical axis by constant amounts +«W. Taking u=u, 
and w =w, with Ay, By, A and B all real and introducing the boundary conditions 
U=U on X= +/ and w=+W on X= +1, we obtain 


U = A, elV3cot?a) Ya a B e—(Vseot?a) la __ Aa tan2a 
0 0 6t 5) 
U = Agen lBeotta) tla 4B eliBeot*a) a ete, 
oF (14.24) 
—— 1 | e(V3cot? a) la aes B e7 (Vacot* a) la Bey a8 B 
a Vee Ve Hie 6t 
= es 1 4,7 W8eot*a) Ia aie ze B, eliootta) Ha __ AG; any 
y3 6T 
These equations yield 
ji en eee 1U + aW tan?« 
e 21 cosh [()/3 cot? «) J/a] — (2/)/3) a tan® sinh [(]/3 cot? «) d/a] ’ 
14.25 
Pes bey 2W cosh [(y3 cot? «) //a] + (2//3) U sinh [(//3 cot? «) /a] 


21 cosh [(j3 cot? «) t/a] — (2/3) a tan? «sinh [(//3 cot? «) //a] ’ 


and B=0. With these values for A,, By, A and B the displacement components 
throughout the membrane are given by (14.21) and (14.23) by taking w =u, and 
wW = Wy, together with v =0. 


Case n=1. When  =1 equation (14.3) becomes 
(a tan «)§ v() + (a tan a)4 (1 — 3 cot?) vf) + (a tan a)? (2 —cot®a) vy =0. (14.26) 


This equation has the general solution 


Ve Ay ef z/atane | Bye" Bia taney, eleciatenc.} Beeq tiene A ai EDR, (14.27) 
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where A,, B,, Az, B,, A and B are integration constants and yw? and ys are 
solutions for A of the equation 


A2 + (1 — 3 cot?a) A+ (2 — cot?) =0. (14.28) 
Then, from (14.27) and the second of equations (14.6), with »=1, we have 


oy tacota [4 (A, els +/a tana nes Bye” salt wes ne 


14.29 
+ be (A, eremia ne a Bie aaieteany | Be 1a27B. ( ) 


From (14.27) and the first of equations (14.6), with 7 =1, we have 


= — La [(ui at 1) (Aggie a Bye" lad =e 


oe (i wis 4) (A, ehstietane Bis pg arene) =f (A ae B2Z)] ‘ (14.30) 

As an example, we consider a problem in which the edges 2 =/ and X= —/ 
undergo equal and opposite small rigid motions. Let the rigid motion of the 
edge S'=/ consist of rotations through angles ew, and ew, about the radial line 
through 9 =0 lying in the plane of the edge and the radial line through O =}2a 
lying in this plane and of translations e/, and eh, respectively in these directions. 
Then, on 2’=/, we have, neglecting terms of higher degree than the first in e, 


u=h, cos 2 + hy sin a5) 
a a 


v= — ss sin 2 da We lice (14.31) 
and a 
w=a (oo ante @ COS 2) { 
a a 
ahus,-on:Li=2; 
u=Al(h,—th,) eo"), 
ey ew i@la 
a v als (ih, +h.) e (14.32) 
w= Rl[—alia,+a,) e%*). 
Similarly, on 2’ = —1, 
u = — Bl (hy — thy) e9%"), 
v= als thy + hq) ef 4 14. 
zis a (le, + ha) oF] (14.33) 
w = — B[— alto, +a,) e%]. 
Introducing 
Uy, =h,— the, vy = —— (thy + hy) 
and w,=—a(t@m,+q@,) on L=/1 
and face (14.34) 
m=—(y—th), %=——(th, +h) 


a 
and w,=a(ia,+o,) on Y=—], 
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we obtain from (14.27), (14.29) and (14.30), with the notation 


Y,=fhl/atana and y,=—pm,l/atana, 


— ta[(ui +1) (4ye"+ Bye”) + 
tei eA Arete Be Mh) teed.i- BU) sh, — a hy: 


— 4a[(ux 4) (Aye*-+ Bem) 4 
+ (42 +1) (Age™+ Bye”) + A — Bl] = —h,tihg, 


A,e™+ Bye ™+ A,e”@+ Boe’+ A+ Bl=— (ih, +h,), 


1 
a 


A,e + Byie+ A,e + Be” | A Bl= - (th, } hs) 


tacot a[u,(A,e%— Bye~™) + ug (Age — Bye-”)] + 


+74? B= — a(t@, + @2) 
and 
tacota[m,(A,e~" — Bye”) + us(Age—” — B, e”)] + 


+1@B=a(io,+q,). 
From (14.36) we obtain 


A,+B,= a tan a(— @, +7@,) cosh yp, , 


My 
Aine Da = raat tan a(w, —7@,) cosh y,, 
Oli ; : 
A,— B,= a (hy — thg) sinh we , 
a uit : ; 
ie eel 8 a 7 (hy — thg) sinh yy, 


= tan « (@, — 7@z2) (wz — uj) cosh y, cosh yp, 
ts fg D (ay 2) (2 — /1) 1 Yo 
eptiana= = “Ay (/ 


ty — thy) (Uy cosh w, sinh ys — fy sinh y, cosh ws) , 


where D and £ are defined by 


4 D = wu, cosh wz sinh py, — “, sinh py, cosh y, , 
an 
E =a(uz — p34) sinh y, sinh y, — 


— [ty Mg! cot «(My cosh y, sinh yp, — 4, sinh y, cosh ws) . 
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(14.35) 


(14.36) 


(14.37) 


(14.38) 


15. Smail Deformation Superposed on Cylindrically Symmetric Large 


Deformation 


The governing equations in this case are equations (11.7), (11.9) and (11.10). 
We shall consider that 7, and #, in these equations, may be expressed as Fourier 
series in @ with period 2za, the coefficients in which are functions of 2’ only. 
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We accordingly write 


co lee) 
1, => fo "o" and, p= 20 (15.1) 
n=0 n=0 
where 7(* and #, are complex functions of ¥ only. As in the case of the small 
deformation of a cylindrical membrane discussed in § 14, if we exclude the 
possibility of dislocations or tears in the membrane, uw, v and w must be continuous 
periodic functions of @ with period 27a. We can therefore express u,v and w 
in the forms 
CO co lee) 
w= nce, p= eo and. w= enor (15.2) 
n=0 n=0 n=0 
where u,, v, and w, are complex functions of 2’ only. We assume also that T, 
and Tt, may be similarly expressed by 


7 =D, eco “and = ee (15.3) 
We therefore consider a solution of equations (11.7), (11.9) and (411.10), in 


which /,, and # are given by 


fe Le ae Se (15.4) 
which has the form 
WU, CV “Geo CPUS a = wg ee 
e (15.5) 
Ty == 70) eras and a — 7?) (EAN 


Introducing (15.5) into (11.7), we obtain 


Y ‘ 5 
‘ Ge (Mn + in7V,) costa + (r’ u,, +: 2' w,) sin? a = 0 
an 


(15.6) 


in(r'u, + 2'w,) +rv, =0. 


Again, introducing (15.4) and (15.5) into equations (11.9), and employing the 
relations (4.6) and the notation 


mW) +=, Hoo, «dodrjad. oc dala, 57) 


we obtain 


in Hes. . 
“(= T,, COS & + 9;, sin a) + 276, sina + 


gee (u, +inry,) cosa + sas ae Al a 0) 0 
ee a COS « i 


yd 1 ota [tn : ; 
(7212 ) (“Fon cosa + x, sim a) + 2(0" 9!" 2" 2") 1 sim a + 


(15.8) 


ae d , / , / , . , , 
= 47 sina qs (ru, + 2'W,) + 47’ sina (7’ uw, + 2’ w}) + 
h(y’*#+272) 7 7.2) 


asin « id Teak 
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Also, introducing (15.4) and (15.5) into (11.11), we obtain 


T 


QZ n2 2 - 9 ” 


a 


2vv' (| n? 2 Rep ee alll 
| 0s a, — sino 
a a 
22’ 
ge om (x sin? @ — - | COS? of at) ty fF 
2in 2 2 
se etd riy’ sin % = =F cos a) v, + 


2Y On if ¢ 
He (r’ sin? «= — cos? - a 
a a 
Quee if 1 n 
SS SN (Zu, a aes ti fe n} 
and (15.9) 


R= 2hp—, [r (rw, +2’ wy) + (+ 2'%) (w, + inr,)] + 


Se et ie 
ats hPn 2 (7 reg *). 


With these expressions for P +Q and R, we have equation (11.10), which, 
the notation (15.7), may be written 


t(P+0)+ \@ 7 — of! 20) sin? a — ai z' costa|t, +R=0. (15.10) 


If 7, z, f, Z and #, are known functions of &, then the five equations (15.6), 
(15.8) and (15.10) can be used to determine w,,, v,, w,,, T, and o,, as functions of 
»' throughout the membrane if appropriate ete aon are imposed on 
the values of these quantities. In this paper, we shall consider, as an example, 
the case when f9)=0 and f”) and #, are real functions of &. Then, from the 
symmetry of the problem, or otherwise, it is easily seen that u,, w, and T,, are 
real, while n and On are imaginary. 

If u,,, %,, Wy, Wy, T, and zo, are known for some value of 2’ (=2, say) then 
they can be ecto pained at the neighboring point X'=2|+A by the following 
procedure: 

(i) from the values of wu, and w), at XY = we can determine the value of w, 
at Y=2,+ Ay from the formula 


[4 Je», p4r= [Uy ]v—», a [un parr Ad; 


(ii) in a similar manner, from the values of w, and w;, at 2 =2,, we can 
determine the value of w, at X=.2,+142; 

(iii) from the first of equations (15.6), we can determine 7v, at 2 =, and 
from the second of equations (15.6), we can determine iv, at Y=}; 

(iv) we then determine 7v, at 2 =2+A2 from the values of 7v, and iv, at 
PEP Fi 

(v) from the first of equations (15.8) we determine io, at X=}; 

(vi) from the values of io, and io, at XY =2), we calculate the value of 1a, 
at Sy = 274-42; 
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(vii) differentiating the first of equations (15.6) with respect to >’, we obtain 


. (tt, + 2in70,) costa + (ru, +2’ w,) sin?a + 
ae 5 (u,, + inrv,) cosa + (ru, +2’ w,) sinra = 0, 
which together with (15.10) provides a pair of equations for the determina- 
tion of uw, and w, at Y=2,; 
(viii) from the values of #, and wu,’ at X =2,, we calculate the value of w,, at 
P=, - AL: 
(ix) from the value of w, and w), at Y=, we calculate w, at Y=2,442; 
(x) from the second of equations (15.8), we calculate 7, at XY = 2); 


(xi) from the values of t, and 1, at 4 =2,, we calculate the value of t,, at 


Se A 


It is seen that if u,,u,,w,,W,,T, and to, are given on the edge S'=/ of 
the membrane, then w,,,7v,,w,, T, and to, can be calculated throughout the 
membrane. If we obtain six linearly independent solutions of this kind deter- 
mined by six different sets of values of u,, u,,,W,, Wy, T, and io, on Y =/, we 
can, by linear combination of these solutions, construct a solution satisfying 
various other types of edge condition, e.g. u,, v,, and w, might be specified on 
Pieviand Dia =<) 


n> 
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